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Computer Assisted 
methods for diagnosing diseases 

Cross Reference to Related Applications 

The present application is a continuation-in-part of co-pending U.S. 
Application Serial Number 08/315,851 filed on September 30, 1994, which 
is a continuation of U.S. Application Serial Number 07/990,772 filed on 
December 14. 1992, now abandoned, which is a continuation-in-part of 
U.S. Application Serial Number 07/964,486 filed on October 21, 1992, 
now abandoned, which is a continuation of U.S. Application Serial Number 
07/806,980 filed on December 12, 1991, now abandoned. 

Technical Field 

The present invention relates to methods for diagnosing, screening 
or prognosing diseases. More particularly, the present invention relates to a 
method for diagnosing, screening or prognosing diseases in himians or 
animals, and for determining the severity and cause of the disease, by 
measuring blood serum levels of specific, predetermined blood constituents 
and then calculating a diagnostic index based on the relationships between 
those blood constituents. Optionally, demographic information, such as the 
age, race or sex of the patient, as well as the patient's medical history can be 
factored into the diagnosis. 

The present invention further relates to an in vitro, serum-based 
computer assisted method for diagnosing, screening or prognosing 
diseases, utilizing a neural network to obtain a diagnostic index. In 
preferred embodiments of the present invention, the method is used to 
diagnose osteopenia and certain cancers, including but not limited to 
ovarian, breast, testicular, colon and prostate cancer. 
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Background of the Invention 

As used herein, the term "disease" is defined as a deviation from the 
normal structure or function of any part, organ or system of the body (or 
any combination thereoO- A specific disease is manifested by characteristic 
symptoms and signs, including both chemical and physical changes. 
Ceruin characteristic signs and symptoms can be quantitated by clinical 
chemical analysis to yield important diagnostic information. For purposes 
of this application, the quantifiable signs, symptoms and/or analytes in 
biological fluids characteristic of a particular disease are defined as 
"biomarkers" for the disease. Current diagnostic methods depend on the 
identification and evaluation of these biomarkers, both individually and as 
they relate to one another. 

Generally, the pathological process involves gradual changes that 
become apparent only when overt change has occurred. In many instances, 
pathological changes involve subtle alterations in multiple biomarkers. It is 
uncommon that a single biomarker will be indicative of the presence or 
absence of a disease. It is the pattern of those biomarkers relative to one 
another and relative to a normal reference range, that is indicative of the 
presence of a disease. 

When individual biomarkers do not show a predictable change and 
the patterns and interrelationships among the biomarkers viewed collectively 
are not clear, the accuracy of a physician's diagnosis is significandy 
reduced. Also, as the number of biomarkers significant for a particular 
disease increases, the number of diagnostic patterns increases, 
correspondingly, pattern recognition and accuracy of a physician's 
diagnosis decreases. 

For example, the following biomarkers collectively show 
characteristic changes in the presence of osteoporosis: calcium, phosphate, 
estradiol (follicular, mid-cycle, luteal, or post-menopausal), progesterone 
(follicular, mid-cycle, luteal, mid-luteal, oral contraceptive, or over 60 
years), alkaline phosphatase, percent liver-ALP, and total intestinal- ALP. 
After measuring these biomarkers, a diagnosing physician would next 
compare the measuring to a normal reference range. While some of the 
biomarkers may fall outside the normal reference range, others may fall 
clearly within the normal reference range. In some circumstances, all of the 
biomarker values may fall within a normal reference range. Presented with 
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such data, a physician may suspect that a patient has undergone some bone 
loss, but will be unable to reach a conclusive and meaningful diagnosis as to 
the presence of the disease osteoporosis. 

The characteristic changes in biomarkers associated with some 
diseases are well documented; however, the quantiutive interpretation of 
each particular biomarker in diagnosing a disease and determining a 
prognosis is not well established. The difficulties inherent in formulating a 
diagnosis from the analysis of a set of laboratory data is best illustrated by 
looking closer at conventional diagnostic methods for a specific disease. A 
discussion of the disease osteoporosis follows. 

The term ''osteopenia" as used herein means any decrease in bone 
mass below the normal. The term "osteoporosis'* as used herein means a 
specific form of generalized osteopenia characterized by a decrease in bone 
density, low bone mass, and microarchitectural deterioration of bone tissue. 

Osteopenia encompasses a group of diseases with diverse etiologies 
typified by reduction in bone mass per unit volume to a level below that 
which is necessary for adequate mechanical support. Osteoporosis is the 
result of the gradual depletion of the inorganic portion of the skeleton and 
can be caused by any number of factors. Primary osteoporosis is an age 
related disorder that is particularly common in women and is characterized 
by decreased bone mass in the absence of other recognizable causes. 
However, osteoporosis occurs in both men and women. . In women it is 
recognized usually at the 5^ or 6th decade, following menopause. In men 
osteoporosis is often recognized around their 6di or 7^ decade. 

The following is a partial list of some of the categories of individuals 
at risk for developing osteoporosis: 

Post-menopausal women 

Cigarette smokers 

Heavy users of alcohol 

Users of a variety of drugs, such as steroids 

Female nmners and ballet dancers 

Male marathoners consuming two few calories 

Bulimics and anorexics 

People with poor diets 

People allergic to dairy products 
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People affected with cancer 

Fair and slim women 

All men and women over the age of 65. 

In addition to being female, the three most significant risk factors are 
poor diet, lack of exercise, and being postmenopausal. Other risk factors 
which are associated with osteoporosis include Caucasian or Oriental 
ancestry, a fair complexions, and a family history of osteoporosis. 

The onset of osteoporosis may be insidious or sudden, following 
trauma. The most common complaint associated with osteoporosis is back 
pain. Eventually, the pain may spread to the pelvis, the thorax, and the 
shoulders. In the spine, the vertebrae can compress, and the back can take 
on a "bent" appearance. Conditions such as kyphosis (humpback) or 
scoliosis may occur. If the spine becomes deformed, other body parts can 
be affected as well. For example, the ribs can be pushed against the pch is. 
or the stomach can be pushed into the pelvis. In addition to spinal 
problems, osteoporosis can also lead to fractures of the hip, wrist, and nhs 
These fractures can occur with only slight trauma and sometimes with no 
trauma at all. Mazess B., et a/., "Bone Density of the Radius, Spine, anJ 
Proximal Femur in Osteoporosis," 7. of Bone cmd Mineral Research, Vvil 
3, pgs. 13-18, (1988); Riggs B. L, et aL, "Involutional Osteoporosis *. 
//£w Engl. /. Med., Vol. 314, pgs. 1676-1686, (1986). The changes 
associated with osteoporosis are gradual so osteoporosis is often not 
detected in its eariy stages. 

From the research perspective, osteoporosis can be classified as 
either primary or secondary to another disease. Primary osteoporosis is 
further classified as juvenile, idiopathic, postmenopausal (Type I) and 
involutional (Type II). It is now understood that accelerated bone loss 
occurs with cessation of menstruation at the time of menopause and in 
women who have amenorrhea as a result of prolactin producing pituitary 
tumor, anorexia nervosa, or intense long-distance running associated with 
undernourishment. These situations are all accompanied by estrogen 
deficiency which is likely to be a major determinant of the accelerated bone 
loss. Bone loss also occurs when esurogen therapy is withdrawn. 

Bone consists primarily of an extracellular mauix containing (by 
weight) approximately 35% organic and 65% inorganic components. The 
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cells of bone represent a minor component of the bone constituents, yet they 
carry out a major portion of the function of the skeletal system. Bone cells 
help maintain serum calcium concentration within a narrow range to regulate 
mineral homeostasis while also being responsible for the continuous 
formation and resorption of the extracellular matrix, allowing response of 
the skeletal system to the mechanical forces resulting firom physical activiQr. 

The major organic component of the extraceUular matrix is collagen. 
This protein has a rigid rodlike structure and is composed of three alpha 
chains held together in helical fashion by covalent and noncovalent forces. 
Multiple collagen molecules form fibrils, and these fibrils in turn are 
arranged in bundles or fibers. It is the bundles or fibers of collagen that can 
be seen in the light microscope as layers or linear arrays. 

Non-collagen components of bone compose a very small portion of 
the organic matrix of the skeleton. This minor fraction consists of proteins, 
glycoproteins, lipids and mucopolysaccharides. Only a few of these 
components have been carefully identiHed and characterized. Two of the 
proteins that have been isolated and studied are osteonectin and osteocalcin. 

Calcium and phosphorus are the main components of the inorganic 
portion of the skeleton. Initially, calcium and phosphorus are deposited as 
amorphous salts but later imdergo rearrangements into a crystalline structure 
that resembles hydroxy apatite (Caio(P04)6(OH)2). Several other ions, 
including Na« K, Mg, and CO3 in varying proportions, may be found in the 
skeletal hydroxyapatite. If there has been fluoride intake, there will also be 
fluoride in the hydroxyapatite. 

Although considerable effort has been expended in studying the 
mechanism of bone mineralization, it is not fully understood. Several 
theories have been suggested to explain the available data. One theory is 
that calcium and phosphorus ions are present in the extracellular fluid in 
amounts exceeding the solubility product of [Ca] x [P]. These ions are 
kept from precipitating by inhibitors of calcification such as pyrophosphate. 
Because osteoblasts contain large amounts of alkaline phosphatase, it has 
been speculated that the activity of this enzyme facilitates mineralization by 
cleaving the phosphate groups, thus altering the Ca:P ratio in the sites of 
calcification. 



PCTAJS95/01379 



6 



In the clinical syndrome of osteoporosis, the reduction in bone mass 
can be attributed to osteopenia because of a dietary deficiency or absorption 
interference of proteins or Vitamin C. It can also be the result of a deficient 
stress stimulus. Osteopenia can also be caused by osteomalacia, a failure of 
proper mineralization of osteoid resulting from bone calciimi or phosphorus 
deficiency or both. It can also be caused by insufficient absorption from the 
intestine due to lack of calcium or a resistance to the action of Vitamin D due 
to failure of its conversion to the biologically active forms, 25 
hydroxychole-calciferol and 1.25 dihydroxychole-calciferol, formed by the 
liver and kidney, respectively. In addition, it may well be caused by an 
abnormal rate of osteolysis due to parathyroid hormone stimulation of 
osteoclastic activity (osteocasts are hematopoitic in origin, arising from the 
migration of monocytes to bone). In reality, most cases of osteoporosis, 
when carefully analyzed, reveal evidence of the causes of osteopenia. 

Chemical analysis of blood may reveal calcium, phosphorus, and 
alkaline phosphatase within the normal range. However, an isoenzyme of 
alkaline phosphatase may be significantly increased. Increased bone 
resorption seen in osteoporotic patients, which occurs as a result of the 
action of osteoclasts, usually involves the dissolution of both minerals and 
organic matrix eventually leading to increased excretion of urinary 
hydroxyproline. Serum estradiol which is secreted almost entirely by the 
ovary is significandy decreased in these patients. This observation is further 
corroborated where it has been demonstrated that exogenous esuogen 
therapy in pecimenopausal women docs delay the onset of post menopausal 
osteopenia. Most experts agree that estrogen appears to decrease bone 
resorption. Weiss N. S., et al., "Decreased Risk of Fractures of the Hip 
and Lower Forearm with Postmenopausal Use of Estrogen," N. EngL J. 
Med,. Vol. 303, pgs. 1195-1198 (1980); Ettinger B., et al, "Long Term 
Esuogen Replacement Therapy Prevents Bone Loss and Fractures,** Ann. 
Intern. Med.. Vol. 102, pgs. 319-324, (1985). 

It is believed that blood levels of calcium are maintained without 
increasing calcium loss from the bone as long as there are normal amounts 
of esuogen and parathyroid hormone present It is reasonably excepted that 
estrogen antagonizes the effect of parathyroid hormone. Estrogen 
deficiency, as seen in peri- and postmenopausal women, results in an 
increase in the sensitivity of bone to parathyroid hormone. This antagonistic 
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relationship eventually leads to an increase in the resorption of the bone and 
contributes to the development of Type I Osteoporosis. Type II 
Osteoporosis is characterized by reduced calcium absorption, which in turn 
results in increased secretion of parathyroid hormone. Type II Osteoporosis 
occurs in older individuals and is associated with wedge fractures of the 
spine and fractures of the hip. A theory recently advanced by Bamhill, et aL 
suggests that evidence of decreased serum estradiol and increased 
lymphocyte alkaline phosphatase, represents an activated immune system in 
osteopenic postmenopausal women. This finding suggests that "imcovered 
estrogen receptors" may induce an inunune reaction which is responsible for 
one form of osteoporosis. Bamhill S.D., et al., ''Osteoporosis: A Possible 
Autoimmune Etiology," Ann. of Clin. Lab. 5ci., Vol. 17, pgs. 255-256, 
(1987). 

In 1987, Bamhill, et aL proposed an autoimmune etiology for some 
forms of osteopenia. In that publication, Bamhill, et aL showed that 
lymphocyte-derived alkaline phosphatase was present in the blood of 90% 
of severely osteopenic women. See Bamhill, etaL^ "Osteoporosis: A 
possible Autoimmune Etiology''^ Ann. of Clin. Lab. ScL, Vol. 17» pgs. 
255-256, (1987)). The concept of a lymphocyte-derived alkaline 
phosphatase is further described by Griffiths, et al. See Griffiths, J., et aL, 
** Separation and Identification of Alkaline Phosphatase Isoenzymes and 
Isoforms in Serum of Healthy Persons by Isoelectric Focusing", Clin 
Chem, Vol. 32, pgs, 2171-2177. (1987). Thus, osteoporosis is well suited 
to diagnosis by interpretation of laboratory data. 

Osteoporosis poses a major health problem in the United States, not 
only for those persons who are already affected but for individuals whose 
diet, life style, and body build make them more likely to develop 
osteoporosis as they age. Postmenopausal osteoporosis is a common 
disorder that results in substantial morbidity and mortality. As many as 25 
million American women suffer from osteoporosis. Osteoporosis is 
responsible for approximately 250,0(X) femoral neck fractures annually in 
the United States, and these fractures are associated with a 20% mortahty 
rate within 6 months. The risk is particularly high among the elderly, who 
also tend to lose bone as a result of the aging process. Wu K., et aL, "Bone 
Resorption Rates in Physiological. Senile, and Postmenopausal 
Osteoporosis." /. Lab Clin. Med., Vol. 69, pg. 810. (1967). 
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It is estimated from a survey of medical clinics that of those 
individuals living to age 85, 32% of the women and 17% of men will 
fracture hips, weakened by osteoporosis. In addition to the pain and 
suffering caused by these fractures, the monetary cost is great, accounting 
for well over 3.8 billion dollars per year in treatment for fractures of 
osteoporosis. Moreover, six to eight months following hip fracture, about 
50% of the osteoporotic patients are in need of assistance with activities of 
daily living and about 25% require nursing home care. Only about 25% of 
the patients fully recover. In view of the costs associated with osteoporosis, 
as measured in dollars and in human suffering, osteoporosis has 
increasingly been perceived as a serious and disabling disease, warranting 
substantial involvement on the part of clinical investigators, governmental 
agencies, and pharmaceutical industries to develop and evaluate potential 
treatments and early detection techniques. 

The NaUonal Institutes of Health 1984 Consensus Conference 
stimulated those interested in bone density with its recommendations for 
"defining persons at risk, and developing safe, effective and low cost 
strategies for fracture protection," See, Office of Medical Applications of 
Research National Institutes of Health, "Osteoporosis: consensus 
conference," JAMA. Vol. 252, pgs. 799-802, (1984). Since 1984, there 
has been a dramatic change in the way physicians view osteopenia. With 
the development and popularization of highly sensitive techniques such as 
DPA and quantitative computed tomography, physicians are now capable of 
measuring the density of the proximal femur and the lumbar vertebrae. (For 
a general review of the methods currentiy available to measure bone density, 
see AviolU "Metabolic Bone Disease and CUnically Related 

Disorders", W.B. Saunders Company (1990)). 

An early decrease in bone mass can be measured by non-invasive 
assessment of the skeleton by four widely available methods, including 
single photon absorptometry. dual photon absorptometry, dual-energy x-ray 
absorptometry. and quantitative computed tomography. A brief discussion 
of these diagnostic methods follows. 

A device called a single-photon absorptometer (SPA) is used to 
measure bone mineral content, primarily in the forearm and wrist. The heel 
bone can also measured using SPA because the heel bone is thought to be a 
predictor of bone loss in the spine. SPA measures primarily cortical bone. 



which is also affected by osteoporosis, though not to the same extent as 
trabecular bone. 

The technique of dual-photon absorptometry (DPA) provides a 
measurement of the total cortical and trabecular mineral content of the hip 
and spine. DPA uses less radiation than conventional x-rays; but, a scan of 
the spine using DPA still exposes the patient to approximately one tenth of 
the radiation that results from a routine chest x-ray. 

Dual-energy x-ray absorptometry (DXA) provides a measurement of 
the amount of bone tissue in the hip and spine. This technique is now used 
routinely because it is faster than DPA. 

Unlike DPA and DXA, quantitative computed tomography, more 
commonly called a CAT scan, can measure the density of either bone or just 
the trabecular portion. CAT scans unfortunately expose patients to higher 
doses of radiation than any of the other techniques. 

Radiographic absorptometry (RA) is a method for non-invasive 
measurement of bone mineral x-rays of the hand. Radiographs, taken with 
a standard x-ray machine, are sent to a central laboratory for computer- 
controlled analysis. 

As osteopenia progresses, both the number and thickness of the 
trabecular imits decrease, causing fragility of bone and an increased risk of 
fractures. The radiographic manifestations of osteoporosis reflect the 
deficiency of the organic matrix and parallel the gross pathologic findings. 
The cortices are thin and the trabeculae fine and sparse; the skeletal 
structures are therefore more radiolucent than normaL The disease process 
eventually affects almost all of the skeletal strucmres. The principal areas of 
demineralization are the spine and pelvis, especially in the femoral neck and 
head. Demineralization is less marked in the skull and extremities. 
Although useful to detect breaks in the bone, ordinary x-rays are not 
sensitive enough to detect osteoporosis until a large amount of bone tissue 
has already been lost, generally from 25% to 40%. By the time 
osteoporosis can be identified by x-ray techniques, the condition is 
advanced. 

Current standard diagnostic techniques, are not effective for early 
detection of osteoporosis. Changes seen in osteoporosis are very gradual, 
and often go undetected in the early stages of the disease. Osteoporosis is 
often not detected in its early stages because bone mass must be decreased 
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by about 30% to 40% before it is apparent using standard x-ray diagnostic 
techniques. Preventing osteoporosis by detecting early bone loss is far 
better than identifying the disease at relatively advanced stages and 
subsequently attempting to prevent its progression. Once major 
deterioration has occurred and gaps exist between the ends of fractured 
trabecular, no current treatment can be expected to restore the lost bone. 
Thus, therapeutic efforts must be directed toward prevention and early 
recognition of the progressive disease so u^eatment can be instituted before 
essentially irreversible structural damage ensues. Cummings S,R., et al., 
"Should Perimenopausal Women Be Screened for Osteoporosis?", Ann. 
Int Med.. Vol. 104, pgs. 745-751, (1986); Courpron P., "Bone Tissue 
Mechanisms Underlying Osteoporosis," Orthop. Clin. North Am,. Vol. 
12, pg. 513. (1981); Frost H, M.. "Mechanical Determinants of Bone 
Modeling," J. MetaboL Bone. Dis. Rel. Res., Vol. 4, pg. 217, (1983). 

One of the problems with the current methods for diagnosing 
osteoporosis is that the procedures do not give any information about the 
underlying cause of the osteoporosis, making it difficult to prescribe an 
appropriate course of treatment for the patient For example, a common 
cause of postmenopausal osteoporosis is an estrogen deficit, which x-ray 
techniques cannot measure. Another problem inherent in the current 
diagnostic methods for osteopenia is that all of the current methods require 
expensive, sophisticated medical instrumentation to perform the bone 
density measurements. Additionally, patients must be exposed to x-rays. 
This makes a general screening of high risk populations impractical due to 
the expense and unavailability of the necessary instrumentation to the 
avOTge clinic. 

In view of the difficulties associated with extracting a diagnosis from 
the laboratory data for a set of predictive biomaikers, there is need for 
automated diagnostic systems that are capable of complex pattern 
recognition. There have been several attempts at using computational 
models to achieve pattern recognition in diagnostics. One of the most 
popular computational method for making diagnoses from multivariate 
laboratory data has been discriminate function analysis. However, 
diagnostic systems that rely exclusively on classical pattern recognition 
technology (geometric, syntactic, template, statistical) are not effective for 
evaluating the characteristic biomarker patterns of many disease states. 
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There is no clear set of roles that accurately describes how to analyze a set of 
biomaiicers to reach a diagnosis. 

In recent years, neural networks have been gaining popularity as a 
means for recognizing and analyzing subtle diagnostic patterns in 
multivariate laboratory data. Neural networks possess the ability to discern 
patterns and trends too subtle or too complex for humans and conventional 
computational methods to identify. While humans can't easily assimilate 
more than two or three variables at once, neural networks can perceive 
correlations among hundreds of variables. Examples of areas in which 
neural networks have been explored for their value in clinical diagnosis 
and/or prognosis include: 

• psychiatry (See Mulsant, B.H., "A Neural Network as an 
Approach to Clinical Diagnosis'*, MD Confuting, Vol. 7. pp. 25-36 
(1990)); 

• autism (See Cohen, et al, "Diagnosing Autism: A Neural 
Net-Based Tool", PCAI, pp. 22-25 (May/June 1994); 

• pediatric radiology (See Boone, JLM„ et al., "Neural Networks 
in Radiologic Diagnosis. I. Introduction and mustration". Invest Radiol, 
Vol. 25. pp. 1012-1016. (1990) and Gross, G.W., et a/., "Neural 
Networks in Radiologic Diagnosis, n. Interpretation of Neonatal Chest 
Radiographs", Invest Radiol, Vol. 25. pp. 1017-1023 (1990)); 

• breast cancer (See Astion, M.L., et al, "Application of Neural 
Networks to the Interpretation of Laboratory Data in Cancer Diagnosis", 
Clin. Chent, Vol. 38. No. 1. pp. 34-38 (1992); Yuzheng, Wl. et al, 
"Artificial Neural Networks in Mammography: Application to Decision 
Making in the Diagnosis of Breast Cancer", Radiology, Vol. 82, pp. 81-87 
(1993); Kappen, H.J., et al, "Neural Network Analysis to Predict 
Treatment Outcome", Annals of Oncology, Vol. 4, Supp. 4, pp. S31-S34 
(1993); and, Ravdin, P.M., et al, " A practical application of neural 
network analysis for predicting outcome of individual breast cancer 
patients". Breast Cancer Research and Treatment, Vol. 22. pp. 285-293 
(1992)); 

• ovarian cancer (See Wilding, P., et al, "Application of 
backpropogation neural networks to diagnosis of breast and ovarian 
cancer". Cancer Letters, Vol. 77, pp. 145-153 (1994)). 
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• thyroid disease (See, Sharpe, et, aL\ "Artifical Neural 
Networks in Diagnosis of Thyroid Function from in Vitro Laboratory 
Tests," Clifu Chem., Vol. 39. No. 11, pps. 2248-2253 (1993)). 

• cervical cancer (See U. S. Patent No. 4,965,725 to Rutenberg); 

and, 

• cardiology (See U.S. Patent No. 5,280,792 to Leong et ai and 
Furlong, J.W., "Neural Network of Serial Cardiac Enzyme Data: A Clinical 
Application of Artifical Machine Intelligence". Clin. Chem., Vol. 96, No. 1, 
pp. 134-141 (July 1991). 

Neural networks are capable of pattern recognition particularly suited 
to making diagnoses. Unlike current methods for arriving at a diagnosis 
from a logical set of rules, neural networks do not require explicit encoding 
of process knowledge in a set of rules. Neural networks learn from 
examples. 

The present invention provides a new method for diagnosing 
diseases and screening patients at risk for a disease, utilizing laboratory data 
to obtain a diagnostic index. It further provides a method for early detection 
and diagnosis of a disease. As applied to osteoporosis, it generates a new 
level of interest in screening individuals who are at risk for osteoporosis 
because this novel diagnostic procedure provides information about the 
underlying cause of the osteopenia. The present invention also provides a 
simple, inexpensive, and rapid method of estimating bone density using 
serum blood levels of specific, predetermined blood constituents. As 
applied to cancer, it provides a simple, inexpensive, accurate and rapid 
method for detecting cancer, resulting in earlier and less invasive treatment. 

Summary of the Invention 

The present invention is directed to an in vitro^ biological fluid-based 
method for diagnosing, screening or prognosing for diseases in which 
laboratory data are manipulated to obtain a diagnostic index. The term 
"biological fluid" includes, but is not limited to. blood, urine, saliva, 
cerebral spinal fluid, synovial fluid, and tears. In the preferred 
embodiments of the present invention, the method of the present invention is 
used to provide a of cancer and osteoporosis. 

The present invention is a method for diagnosing, screening or 
prognosing a disease in a human or animal comprising the steps of 
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measuring the concentrations of a predetermined set of biomarkers known to 
be associated with the disease from a biological fluid from the human or 
animal; scaling the digitized values of the analytes; and sending the scaled 
values to a trained neural networic, whereby the output values from the 
neural network tend toward the upper value when the human or animal has 
the disease and the output values tend toward the lower value when the 
human or animal does not have the disease. 

The present invention also comprises an apparatus for diagnosing, 
screening or prognosing a disease in a human or animal comprising a means 
for digitizing the concentrations of a predetermined set of biomarkers known 
to be associated with the disease from a biological fluid from the human or 
animal; a means for scaling the digitized values; and a trained neural network 
coupled to the digitizing and scaling means for generating network output 
values between an upper and lower value, whereby the output values tend 
toward the upper value when the human or animal has the disease and the 
output values tend toward the lower value when the human or animal does 
not have the disease. 

The present invention also comprises a method for diagnosing, 
screening or prognosing a disease in a human or animal comprising the 
steps of measuring the concentrations of a predetermined set of biomarkers 
known to be associated with the disease from a biological fluid from the 
human or animal, scaling the digitized values of the analytes, and 
introducing the scaled values to a first trained neural network, whereby the 
output values from the first neural network tend toward the upper value 
when the himian or animal has the disease and the output values tend toward 
the lower value when the human or animal does not have the disease; and 
sending the output value from the first neural network and a second set of 
predetermined biomarkers, which could include one or more of the 
biomarkers in the first set of predetermined biomarkers, to a second trained 
neural network, whereby the output values from the second neural network 
tend toward the upper value when the human or animal has the disease and 
the output values from the second neural network tend toward the lower 
value when the human or animal does not have the disease. 

In accordance with the first embodiment of the present invention, a 
trained neural network is utilized to determine a diagnostic index 
corresponding to the presence and severity of a disease by analyzing a set of 
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predetennined biomarkers for that disease. In accordance with the 
invention, the concentrations, or in some cases* the presence of certain 
biomarkers related to the incidence of a particular disease are determined for 
a patient These values are then scaled and the scaled values are then sent to 
a trained neural networic to yield a diagnostic index. A neural network is 
trained by introducing a population of patients in which a disease state is 
known, along with the biomaiker values for those patients and ''teaching*' 
the neural network to recognize the patterns in the biomarkers. After the 
neural network is trained, biomarker values from patients with unknown 
disease states are introduced to the trained neural network. The neural 
network then processes the information to produce a value corresponding to 
a diagnosis of the presence or absence and the severity of a particular 
disease. 

The inventors propose that the artificial neural network, especially 
the multi-layer feedforward network, may, through their weighi 
connections, correspond to biomarker patterns that are important for 
categorizing diseases. Additionally, the inventors propose that the ncur^il 
network is able to identify unique patterns of biomarkers associated with a 
variety of disorders that may help to classify borderline cases that do noi 
appear to fit into either a malignant or benign pattern. 

A second embodiment of the present invention involves a two step 
analysis of the biomarkers by neural network. This avoids the bias created 
by a dominant predictive variable when training a networic. The dominant 
biomarker or predictive variable is excluded from the first analysis by neural 
network and then included in a second analysis by neural network. For 
example, if age is thought to be the dominant predictive variable in the 
diagnosis of osteoporosis, that variable is not included in the training of the 
flrst neural networic, and the training data set is limited to the other selected 
biomarkers. After obtaining a diagnostic index using the first set of 
biomarkers, a second neural network is trained using the diagnostic index 
and the entire set of input variables, including age, to yield another 
diagnostic index. The final diagnostic index is a composition of an artificial 
neural network generated index and results from heuristic analysis using 
other non-numerical patient information. 

The present invention further comprises an apparatus capable of 
diagnosing, screening or prognosing a disease comprising a sample 
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receiving means, a sample detecting means capable of detecting the quantity 
of one or more biomarkers in a biological fluid, and, after digitalizing the 
quantity of analyte in the biological fluid, scaling the digitized values and 
introducing the scaled data into a neural network which yields a therapeutic 
index on a printer or, optionally, on a video display. 

A third embodiment of the present invention includes a method for 
diagnosing and determining the severity and underlying cause of osteopenia 
using blood concentrations of a specific, predetermined set of blood 
constituents. The method comprises the steps of (1) measuring the severity 
of the disease in a set of humans or animals with varying severity of disease 
by a standard method, (2) assigning the severity of disease a numerical 
value on a severity scale, the scale being from no disease to severe disease. 
(3) measuring the blood concentrations of the predetermined set of blood 
constituents in the set of humans or animals with varying severity of the 
disease, and then (4) calculating a numerical relationship between the set of 
blood concentrations and the severity of disease. 

The present invention also comprises a simple and rapid method of 
determining severity of osteopenia in a patient In a preferred embodiment 
the method comprises determining the serum level of the following scrum 
constituents: calcium, phosphate, total alkaline phosphatase, an alkahnc 
phosphatase isoenzyme, estradiol, and progesterone. The alkahnc 
phosphatase isoenzyme is preferably t-lymphocyte derived alkaline 
phosphatase or blood, liver or intestinal alkaline phosphatase isoenzyme. 
The results of these tests are then introduced into an algorithm. Optionally, 
the age of the patient may also be factored into the equation. The bone 
density coefficient that is calculated by the algorithm correlates to a very 
high degree to bone density as measured by standard methods, such as 
radiographic absorptometry, quantitative computed tomography, dual 
photon absorptometry and diiea measurement of bone doisity. The bone 
density coefficient that is measured is then compared to an osteopenic 
severity scale. 

Using the six seium constituent concentrations, the present invention 
can be used to determine the osteopenic state of a patient as well as give an 
indication of the underlying cause of the osteopenia. The present invention 
can be correlated to any method of measuring bone density simply by 
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recalculating the coefficients in the algorithm using a multiple linear 
regression analysis. 

Accordingly, it is an object of the present invention is to provide a 
method for diagnosing, screening or prognosing and determining the 
severity of a disease using the concentrations and values of a predetermined 
set of biomaikers. 

It is a further object of the present invention to provide a more 
accurate method for screening, prognosing and diagnosing breast cancer 
using biological fluid-based laboratory data. 

It is another object of the present invention to provide a more 
accurate method for screening and diagnosing ovarian cancer using 
biological fluid-based laboratory data. 

Yet another object of the present invention is to provide a more 
accurate method for screening, prognosing and diagnosing colon cancer 
using biological fluid-based laboratory data. 

Furthermore, it is an object of the present invention to provide a 
more accurate method for screening, prognosing and diagnosing prostate 
cancer using biological fluid-based laboratory data. 

Additionally, it is an object of the present invention to provide a 
more accurate method for screening, prognosing and diagnosing testicular 
cancer using biological fluid-based laboratory data. 

Another object of the present invention is to provide a method for 
diagnosing cancer which will provide a better understanding of the probable 
cause of the cancer. 

Another object of the present invention is to provide a diagnostic test 
for cancer which can be used to screen large numbers of individuals. 

Another object of the present invention is to provide a simple and 
rapid chemical test for diagnosing osteoporosis. 

Another object of the present invention is to provide a test for 
osteoporosis which will also give information as to the underlying cause of 
the osteopenic condition. 

Another object of the present invention is to provide a diagnostic test 
for osteoporosis which can be used to screen large numbers of individuals. 

Yet another object of the present invention is to provide a method for 
diagnosing osteoporosis and determining the underlying cause of the 
osteopenia without having to subject the patient to radiation. 
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These and other objects, features and advantages of the present 
invention will become apparent after a review of the following detailed 
description of the disclosed embodiment and the appended claims. 

Brief Description of the Figures 

Fig. 1 illustrates a feed forward neural network having multiple 
outputs. 

Fig« 2 illustrates a feed forward neural network having a single 

output. 

Fig. 3 is an equation illustrating the mathematical relationship 
between the input and output of a neuron. 

Fig. 4 is a schematic illustration of the second preferred 
embodiment of the present invention. 

Fig. 5 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the present invention over time 
for a ovarian cancer patirent 

Fig. 6 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the present invention over time 
for a prostate cancer patirent 

Fig. 7 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the present invention over time 
for a colon cancer patirent 

Fig. 8 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the presmt invention over time 
for a breast cancer patient 

Fig. 9 shows the patient data used to develop and test the 
osteoporosis neural network diagnostic system. 

Fig. 10 shows the patient data used to test the osteoporosis neural 
netwoik diagnostic system. 

Fig. 11 is a program list of the algorithm of the third embodiment 
of the present invention. 

Fig. 12 shows the discrepancies between the DPA measurements 
and the measurements obtained using the method of the present invention 
for individual subjects. 

Fig. 13 shows the correlation of the results obtained using the 
present invention and DPA. 
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Fig. 14 shows the normality of values predicted by the present 
invention. 

Fig. 15 shows the normality of residuals (differences between DPA 
measurements and values predicted according to the method of the present 
invention), as evidenced by the arrangement of points along a diagonal 
straight line. 

Fig. 16 shows that the variance of residuals. 

Fig- 17 shows the independence of residuals. 

Fig. 18 shows that virtually all observed DPA measurements are 
predicted by the method of the present invention. 

Fig. 19 shows the training data used to construct the prostate 
cancer neural network prognostic system. 

Fig. 20 shows the relationships between the input neurons, hidden 
layer and output neuron of a neural network in one embodiment of the 
present inventioru 

Detailed Description 

The first embodiment of the present invention is directed to an in 
vitro, serum-based computer assisted method for screening, prognosing and 
diagnosing diseases utilizing a neural network to obtain a conclusive 
diagnosis. The present invention can be adapted to existing diagnostic 
devices that have a collection means, a sample detecting means capable of 
detecting the quantiQr of an analyte in a biological fluid and a means of either 
printing or displaying the results of the tests on video display means. 

The inventors have discovered that biomarkers collectively alter in 
response to a disease process, and collectively constitute a new diagnostic 
biomarker with better disease predictability than the individual biomarkers. 
When the biomarkers are processed and analyzed as a group to yield a single 
diagnostic index, the sensitivity and specificity of the diagnosis is increased, 
making it possible for a physician to detect the presence of a disease earlier 
and with greater precision, or estimate a prognosis with greater precision, 
than by analysis of the individual biomarkers. 

In accordance with the method of this invention, blood is first 
collected from a patient, and then centrifuged to separate the red blood cells 
from the serum. It should be understood that biological fluids other than 
blood can be used in practicing the present invention. Next, the serum is 
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analyzed, using standard laboratory techniques* to determine the 
concentrations, or in some cases the presence or absence, of specific 
predetermined biomarkers related to the incidence of a particular disease. It 
is to be understood that this process can be carried out automatically in 
conventional diagnostic machines. For purposes of illustration, a 
description of the methods for obtaining the values for the biomarkers for 
osteopenia is provided elsewhere in this sectioiL 

The biomarkers relied upon to diagnose a disease by the method of 
the present invention must be predictive of the suspected disease and must 
be statistically significant for analysis by a neural network. The selecdon of 
biomarkers that offers statistically significant discriminating power in the 
diagnosis of disease involves several steps. First an inventory of 
biomarkers that have shown certain relevancy in the diagnosis of the disease 
of interest must be conducted. In general, only the biomarkers that reflect 
different aspects of the disease process or other diagnostic information need 
to be included. Second, the selected biomarkers need to have a reasonable 
diagnostic value in terms of sensitivity, specificity, and positive and 
negative predictive powers. The design and implementation of experimental 
protocol from which the biomarkers are developed and evaluated should 
also be considered. Third, if the number of candidate biomarkers is large, a 
formal discriminating power analysis may be conducted. However, many 
of the standard statistical analysis methods may not be adequate for highly 
nonlinear classification problems. If the number of candidates are not too 
large, they may be all included in the initial attempt of neural network 
training. If one or several of the input biomarkers to the network are 
irrelevant to the classification decision making process, it will be reflected in 
the network connection weights of the trained neural networks. These 
values may then be removed from the biomarker set for a particular disease. 
Other methods for evaluating the statistical signiHcance of a biomarker 
selected for analysis by neural network and selecting biomarkers for training 
a neural network are well known in the art 

Biomarkers which meet the criteria delineated above, namely, they 
are predictive of a particular disease and statistically significant for analysis 
by neural network, are identified below for a sample of diseases including 
ovarian cancer, prostate cancer, colon carK:er, breast cancer, testicular cancer 
and osteoporosis. 
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Ovarian Cancer 


Prostate Cancer I 


Colon Cancer 


La^a-p® 


LASA-P® 


LASA-P® 


CAiis 


PAP 


CA 19-9 


dWjOk 


pSa 


Cea 



Pr ostate Cancer II 

pAI^ "" 
MA 

Acid Phosphatase " 



Breast Cancer 


Testicular Cancer 


Osteoporosis 


LASA-P(g> 


LASA-P® 


Calcium 


CfeA 


AFP 


Phosphate 


H£R2/neu 
in Plasma 


HCG'Beta 


Estradiol 




Ca i5.3<8> 


Pro^sterone 






ALP 






ALP Isoenzyme 1 






ALP Isoenzyme 2 



A key to the abbreviations used above is provided below: 



AFP: 


Alpha-Feioproiein 


CA125: 


Canco" Antigen 125 


Ca 15-5*** 


Bieast Antigms 1 15D8A3F3 


Ca 1^9: 


Carbohydrate Antigen 19-9 




Caicinoembryonic Antigen 


CK-BB: 


Creatinine kinase, BB subfraction 


DMHOk: 


Ovarian marker NB/70K 


HCG-Beta: 


Human Chorionic Gonadotropin, Beta Sub- 
Unit 


HER 2/neu in 

Plasma: 


c-erb B-2 (H£R2/neu) oncoprotein in plasma 


LASA-P®*: 


Lipid- Associated Sialic Acid in Plasma 


PAP: 


Prostatic Acid Phosphatase 


PSA: 




fic Antigen 



*LAS A-P is a registered trademark of DIANON Systems. Inc. 
10 **CA 15-3 is a registered trademark of Centocor. Inc. 
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A wide number of diseases may be diagnosed in accordance with the 
method of the present invention. To be suitable for diagnosis by the present 
method, biomarkers for the disease must be quantifiable. The biomarkers 
must also be predictive of the disease and must be statistically significant 
relative to one another. The method of the present invention is equally 
suited to the diagnosis of any disease in which serum biomarkers can be 
identified, including but not limited to infectious diseases, and genetic 
abnormalities. 

After determining the serum biomarkers for a disease, the biomarker 
values are analyzed by a trained neural network to yield a single diagnostic 
value. The most common neural network architecture for pattern 
classification problems is the feedforward network, which typically consists 
of an input layer, one or more hidden layers, and an output layer. Figs. 1 
and 2 illustrate the arrangement of neurons in two different feedforward 
networks. 

The elements that make up each layer of a neural network are 
refened to as neurons or nodes. Inputs are fed forward from the input layer 
to the hidden layers and then to the output layer. The number of neurons in 
each layer is determined before the networic is trained. Typically, there is 
one input neuron or node for each input variable, and one output node for 
each output The inputs to the neural network are predictor variables. 
These predictor variables can be quantitative or qualitative. Neural networks 
make no data distribution assumptions and can simultaneously use both 
quantitative and qualitative inputs. In the present invention, the biomarker 
values, after a linear prescaling to values between 0.0 and 0.1, constitute the 
input variables. 

The outputs of the network represent output categories. For 
example, a malignancy may be represented by maximal output of the 
malignant output neuron and silence of the boiign neuron, whereas a benign 
process is represented by maximal output of the benign neuron and silence 
of the malignant neuron. A simple arithmetic function combines outputs of 
the two neurons to yield a single diagnostic index. In the alternative, a 
single output netu-on may be used. An output of greater than 0.5 would 
indicate a malignancy and an ou4>ut of less than 0.5 would indicate a benign 
condition. In this way a diagnostic index is directly obtained. 
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The number of hidden layers and the number of nodes in the hidden 
layers are configurable parameters that have a significant influence on the 
performance of the network. In practice, the optimal number of hidden 
neurons is determined empirically. The means for deteraiining the optimum 

5 numbers of hidden neurons is well known to those skilled in the art and 

depends on the complexity of the problem being solved 

In the present invention, the main network model is a multi-layer 
feedforward perception using a backpropogation training algorithm. The 
number of hidden layers and the number of neurons in each hidden layer 

10 was determined to adequately match the level of complexity of the diagnostic 

problem. With the assumption that the samples in the training set are 
representative of all possible situations encountered in real applicationswith 
no significant contradictions, and and the number and stratification of 
samples in the generalization test are statistically adequeaie, the criteria 

15 ouUined below is used to determine if a chosen network configuration is 

appropriate. 

If the network continues to fail to correcUy classify large portions of 
the samples in the training set, even after many adjustments of training 
algorithm parameters, the network complexity should be increased. 
20 On the other hand, if the the network acheives a high rate of 

correcUy classifying the training set but fails to accurately classify a large 
number of samples in the testing set, network smicture is probably too 
complex for the problem being solved, Le. it has sufficient inherent 
flexibility to fit the training data set. but not sufficient predictive power to 
25 classify the lest data set If this is the case, the number of neurons in the 

hidden layers should gradually be reduced, or, if there are mutiple hidden 
layers, the hidden layers should be gradually reduced. 

Note that it is not always necessary to have a large uaining sample 
set If the samples in a training set have already represented all possible 
30 cases with adequate sutistical significance, the addition of new samples 

genewrally does not increase the amount of information in the training 
samples. Instead it may decrease the useful information to nosie ratio in the 
samples. At the other extreme, too small a training set will generaUy not be 
able to cover all possible variations in the population. The resultant network 
35 often simply memorizes all the cases in the training set and does not 

generalize at all. 
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The input and output layers are not directly connected. Every input 
neuron is connected to every neuron in the following hidden layer and 
neuron in ahidden layer is connected to every neuron in the following 
adjacent hidden layer or output layer, depending on the number of hidden 
layers. Each of the multiple connections to a particular neuron is weighted. 
In the hidden and output layers, each node sums the input activations, 
multiplied by the respective connection weights plus a bias term. The 
weighted sum then passes through a non-linear output function, typically a 
sigmoidal function, which gives the network the ability to represent complex 
non-linear relationships. A neuron fires if the sura of the weighted inputs to 
it are greater than a threshold value. As illustrated in Fig. 3, once a neuron 
is above a threshold, the magnitude of its output is a sigmoid function of the 
net input The end residt of activity in the neural network is the net ou^ut, a 
complex nonlinear function of the inputs. 

In summary and in accordance with the present invention, first the 
values of the biomarkers for a specific disease are determined and scaled. 
The biomarkers are fed forward from the input layer to the hidden layer (or 
layers) and then to the output layer of the neural network. The number of 
neurons in the input layer is determined before the network is trained and 
correspond to the number of biomarkers predictive for a specific disease. 
The biomarkers are preselected. There is one input neuron for each 
diagnostic variable or biomaiker, and one output neuron for each desired 
output. Other than the identified biomarkers, diagnostic variables may 
include demographic information such as the age, race or sex of the patient 
and whether the patient is pre- or post-menopausal. The number of neurons 
in the output layer depends on the type of output desired. The number of 
neurons in the hidden layer is determined eupencally during training. 

The neural network used for diagnosing a specific disease must be 
trained to do so. In accordance with the present invention, the neural 
network is trained by back propagation. Back propagation refers to the 
technique of training a neural network to accurately model a set of input and 
output examples by determining the best connection weights between the 
values, and is well known in the art Other techniques which may be used 
to train a neural network for purposes of this invention may include any 
other non-linear global optimization technique, such as the genetic search 
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algorithm; however, the feed forward, back propogation network is most 
popular. 

At the beginning of training, the connection weights in the network 
are randomly initialized. The training data is then presented to the network 
one at a time. In accordance with the present invention, the training data 
consists of the biomarker values for a group of patients, and the diagnosis 
for each of those patients. The biomarker values are the input variables used 
to train the network. For each patient, the network uses the patient's 
biomarker values to estimate a diagnosis, which is then compared to the 
actual diagnosis. If the network's diagnosis is correct, then the connection 
strengths and thresholds within the network are not changed, and the next 
patient is presented to the netwoik. If the estimate of the diagnosis is not 
coirect, the connection weights and thresholds in both the hidden layer and 
the output layer are adjusted to reduce the size of the classification error. 
After adjustments are made, the next patient is presented. Training proceeds 
until all patients in the training group are correctly classified. If the network 
is not able to completely classify all the data it will train indefinitely, in 
which case training is tenninated. 

When training the neural network, the trainer may set the 
decision limits regarding the definition of a classification error, i.e. an 
incorrect diagnosis. The relevant parameter is the error tolerance, which 
specifies how close the estimated output has to be to the actual output to be 
correct For example, if two output neurons are used and the training 
tolerance is set at 5%. the estimate of malignancy is considered correct if the 
malignant output neuron fires at 95% of maximum and the benign neuron 
fires at 5% of maximum. Similarly, a correct estimate of a benign diagnosis 
means that the benign output neuron firs at 95% of maximum, while the 
malignant neuron fires at 5% of maximum. The methods for determining a 
classification error are well known to those skilled in the art. 

In the preferred embodiment of this invention, if a single output 
neuron is used, a benign normal diagnosis is set at an output of 0.1 and a 
malignant or abnormal diagnosis is set at an output of 0.9. Error tolerance 
is an adjustable parameter and is significant in determining the success of the 
network at making an accurate diagnosis. It is preferred to use either one or 
two output neurons. 
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After the neural network is trained for the desired disease, biomaiker 
values from patients with unknown disease conditions are introduced to the 
trained neural network. The neural network then processes the information 
to produce a value corresponding to a diagnosis of the presence or absence 
of the particular disease. In accordance with the present invention, this is 
accomplished by using either one single output neuron or two output 
neurons. If two output neurons are used, the ouq)ut from the two neurons 
are combined to generate a single diagnostic index. 

As illustarted by Fig. 4, in a second embodiment of the present 
invention, the diagnostic value obtained by analysis of the biomaikers by a 
trained neural network is further analyzed by a set of heuristic rules in 
combination with additional patient information. The additional patient 
information includes things such as family medical history and demographic 
information. This data is then processed to yield a second single diagnostic 
value. 

This invention is further illustrated by the following examples, 
which are not to be construed in any way as imposing limitations upon the 
scope thereof. On the contrary, it is to be clearly understood that reson ma> 
be had to various other embodiments, modifications, and equivalents 
thereof, which, after reading the description herein, may suggest themselves 
to those skilled in the art without departing from the spirit of the present 
invention and/or the scope of the appended claims. 

In the following examples which utilize a neural network in the 
analysis of the data, a Neuialshell 2, Release 1.5 (Ward Systems Group, 
Inc.) neural network development program was used for the training of the 
neural network on a Pentium 60 mhz computer (Magitronic, Inc.). 

Example 1 

The following example describes the training of a neural network to 
diagnose colon cancer. 

A total of 167 patients ( 48 with diagnosed colon cancer, and 1 19 
normal) were divided into 2 groups, a training set and a generalization 
testing set The training set contained approximately 60% of the patients (90 
patients, 27 with diagnosed colon cance, and 63 normal), and the 
generalization testing set contained approximately 40% of the patients (77 
patients, 23 with diagnosed colon cancer, and 54 normal). 
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The initial network architecture was selected based on the level of 
complextity of the classification task. A multi-layer feedforward network 
was used. Selection of the initial architecture involved the selection of the 
number of hidden layers and the number of neurons in each hidden layer. 
5 Several trial iterations were performed to determine and adequeate 

configuration that showed good results on both the training sample set and 
the generalization test sample set This particular neural network had 5 input 
neurons, one hidden layer having 10 neurons and one output neuron. 

Initially, connection weights among the neurons were randomly set 
10 The neural network had five input neurons, corresponding to five input 

variables for diagnosing colon cancer: LASA-P®; CA 19-9, CEA. Age, 
and Sex. During gaining, the five input variables for each patient were first 
linearly scaled into the continues range between 0.0 and 1.0. The resultant 
five numbers were then presented as an input vector to the input neurons of 
15 the artificial neural networit. 

For each of the input vectors, the network generated an output based 
on the connection weights among the network neurons. The output can be a 
single value or a vector of numbers, depending on the number of output 
neurons used. Each neuron in the network participated in the output 
20 calculation by passing the sum of all inputs to the neuron through a non- 

linear s-shaped function (often a HgisUc function) and sending the result to 
each and everyone of the neurons in the following adjacent layer. The 
generated output was compared to the desired "target" output A value of 
0,1 corresponded to a diagnosis of normal and an output of 0.9 
25 corresponded to a diagnosis of of abnormal. The difference was used to 

calculate an error term to guide the training algorithm, i.e.. the 
backpropogation algorithm, in the adjustment of network connection 
weights in an attempt to reduce the differences between networic outputs and 
target values over the training sample set 
3Q To counteract over learning, which can occur when a neural network 

is too large for the problem and trains too long, the network was tested on 
the generalization test set every 200 iterations. The networks for which 
generalization testing showed reduced root mean squared error over 
previous tests were saved and used in later calculations or generalizations. 
35 This procedure was repeated until generalization failed to improve over 2000 

successive tests. 
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To accelerate the rate of solution during the initial stage of training, 
the learning rate for the hidden layer was initially set at 0.1 and momennun 
was set at 0.05, for the first several iterations. After the first several 
iterations, the learning rate was lowered to 0.01. and momentum to O.OOOS 
to narrow in on the appropriate solution. 

After training, the neural network correctly identified 95% of the 
patients without cancer in the training data set, and 92% of the patients with 
cancer. 

When presented with the generalization test results, the trained 
neural network correctly identified 100% of the people with normal 
diagnosis and 95% of the people with cancer diagnosis. 
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Exampk 2 

The following example illustrates the value of the present invention 
for diagnosing ovarian cancer. 

The bioraarker values for LASA-P®, CA 125. and DM/70K were 
measured in a patient suspected of having ovarian cancer. Two of the 
biomarker values fell within the normal reference range, and one value fell 
outside the normal reference range. The values were scaled and analyzed by 
a trained neural network, yielding a diagnostic index of 75, indicating the 
presence of ovarian cancer. The diagnotic index generated by the neural 
network (between O.I and 0.9) is multiplied by 100 to obtain a whole 
number. The normal reference range for the diagnostic index was 0-50. 
The biomarker values measured, together with the normal reference range 
for those values are provided below: 



Biomarker 


Patient Value 


Normal Reference 
Ranee 


LASA-P® 


1§ 




CAllS 




0-35 U/ml 


DM/70K 


45 


0-35 U/ml 


Neural Net 
Diagnostic bidex 


75 


0-50 



30 



The patient underwent surgery to treat, the cancer. Fig. 5 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 
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Kicample 3 

The following example illustrates the value of the present invention 
for diagnosing prostate cancer. 

The biomarker values for LAS A-P®, PAP, and PSA were measured 
in a patient suspected of having prostate cancer All three of the biomarker 
values feU within the normal reference range. The values were scaled and 
analyzed by a trained neural network, yielding a diagnostic index of 90, 
indicating the presence of prosute cancer. The normal reference range for 
the diagnostic index was 0-50. The biomarker values measured, together 
with the normal reference range for those values are provided below 



Biomaiker 


Patient Value 


Nonnal Refoence Range 


LASA-P® 


16 


15-20 me/dl 


PAP 


1.5 


0-2.5 ng/ml 


PSA 


1.0 


0-2.5 n^ml 


Nexiral Net 


90 


0-50 



The patient underwent surgery to treat the cancer. Fig. 6 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 



Exampk 4 

The following example illustrates the value of the present invention 
for diagnosing colon cancer. 

The biomarker values for LASA-P®, CA 19-9, and CEA were 
measured in a patient suspected of having colon cancer. All three of the 
biomarker values feU within the normal reference range. The values were 
scaled and analyzed by a trained neural network, yielding a diagnostic index 
of 93. indicating the presence of colon cancer. The normal reference range 
for the diagnostic index was 0-50. The biomarker values measured, 
together with the normal reference range for those values are provided 
below: 
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Biomarker 


Patient Value 


Normal Reference Range 


LASA-P® 


17.9 


15-20 rap/dl 


CA19-9 


5.0 


d-37 tjyral 


CEA 




0-2.5 ne/ml 


Neural Net 
Diagnostic Index 


95 


0-50 



The patient underwent surgery to treat the cancer. Hg. 7 graphically 
illustrates the change in the neural network diagnostic index over time» 
before and after surgery. 



10 



15 



Example S 

The following example illustrates the value of the present invention 
for diagnosing breast cancer. 

The biomarker values for LASA-P®, CEA. HER2/neu in Plasma 
and CA 15-3® were measured in a patient suspected of having breast 
cancer. All four of the biomarker values fell within the normal reference 
range, although two of the biomarker values were on the high end of the 
normal range. The values were scaled and analyzed by a trained neural 
network, yielding a diagnostic index of 91, indicating the presence of breast 
cancer. The normal reference range for the diagnostic index was 0-50. The 
biomarker values measured, together with the normal reference range for 
those values are provided below: 



Biomaiker 


Patient Value 


Normal Reference Range 


LASA-P® 


18 


15-20 me/dl 


CEA 


1.0 


0-2.5 ng/ml 


HER2/neu in Plasma 


19 


0-20 U/ml 


CA 15-3® 


8 


0-30 U/ml 


Neural Net 
Diagnostic Index 


91 


0-50 



20 



The patient underwent surgery to treat the cancer. Fig. 8 graphicaUy 
illustrates the change in the neural network diagnostic index overtime, 
before and after surgery. 
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Example 6 

This example illustrates the construction and training of a neural 
network for diagnosis of ovarian cancer. 

Table I below provides the data used to train the network (age, 
LASA-P®. CA 125. and DM/70K). the actual histological diagnosis, and 
diagnostic index determined by the neural networic A histological index of 
0.1 corresponds to a normal diagnosis and a histological index of 0.9 
corresponds to an abnormal diagnosis. A neural net diagnostic index less 
than .5 corresponds to a normal diagnosis. A neural net diagnostic index 
greater than .5 corresponds to an abnormal diagnosis. 



TABLE 1 

Training Data for Ovarian Cancer Diagnostic System 



Patient 
Number 


Age 


^SA-P<S 


CA 125 


DM/70K 


Histological 1 
Diagnosis 


Neural Net 
Diagnostic 
Index 


94 


34 


15.8 


4 


7 


0.1 


0.225 


99 


27 


19.1 


32 


7 


0.1 


0.357 


100 


33 


27 


19 


27 


0.1 


0.317 


101 


32 


16.8 


7 


0 


0.1 


0.236 


102 


29 


17.5 


8 


13 


0.1 


0.186 


106 


25 


18 


3 


0 


0.1 1 


0.137 


110 


36 


17.7 


21 


0 


0.1 


0.419 


112 


37 


12.4 


10 


0 


0.1 1 


0.318 


114 


27 


16.2 


7 


11 


0.1 


0.159 


2001 


44 


43.5 


1880 


182 


0.9 


0.657 


2007 


51 


27 


142 


86 


0.9 


0.75 


2009 


64 


20.4 


55 


32 


0.9 


0.87 


2010 


72 


27.1 


7 


43 


0.9 


1 0.789 


2015 


39 


22 


60 


55 


0.9 


1 0.568 


2019 


35 


16.5 


13 


0 


0.9 


0.327 


2020 


66 


23.1 


19 


46 


0.9 


0.79 


2022 


69 


23.2 


189 


28 


0.9 


0.905 
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91 








OX 


9^ 1 


\ASLf\ 
l*foU 


35 




U.<$4o 


2U32 


TO 


/ 


TO 

to 


21 


0.9 


0.946 


2034 


d5 


3U.O 


43y 


A < 

43 


o.y 


0.876 


2036 


48 


14.2 


13 


0 


0.9 


0.538 


2037 


72 


16.5 


15 


0 


0.9 


0.858 


2038 


68 


13.8 


25 


11 


0.9 


0.851 


2039 


69 


12.1 


10 


3 


0.9 


0,796 


2040 


79 


25.2 


50 


45 


0.9 


0.939 


46 


42 


22.6 


3 


15 


0.1 


0.34 


49 


25 


16.1 


6 


7 


0.1 


0.14 


51 


28 


15.2 


17 


0 


0.1 


0.251 


52 


40 


19.1 


0 


0 


0.1 


0.303 


53 


27 


25 


4 


5 


0.1 


0.17 


56 


31 


24,6 


0 


0 


0.1 


0.196 


57 


33 


13.7 


7 


3 


0.1 


0.235 


58 


25 


24.8 


12 


9 


0.1 


0.19 


59 


25 


25 


9 


11 


0.1 


0.169 


60 


25 


17.4 


5 


0 


0. 1 


0.146 


61 


35 


23.3 


0 


15 


0.1 


0.217 








19 


n 

V 


n 1 


n 990 


63 


28 


18.3 


2 


1 


0.1 


0.16 


65 


33 


15.4 


8 


9 


0.1 


0.236 


68 


30 


22.8 


2 


15 


0.1 


0.171 


70 


28 


15.4 


24 


1 


0.1 


0.305 


71 


31 


20 


6 


19 


0.1 


0.193 


73 


30 


18 


9 


0 


0.1 


0.228 


75 


35 


23.1 


1 


0 


0.1 


0.25 


76 


35 


20.4 


3 


0 


0.1 


0.258 


80 


34 


18.9 


1 


6 


0.1 


0.215 


82 


32 


18.2 


72 


152 


0.1 


1 0.435 




wo 96/12187 PCTAJS95/01379 

32 



84 


33 


17.4 


19 


21 


0.1 1 0.305 


85 


26 


15.6 


7 


5 


0.1 1 0.156 


89 


26 


20.7 


12 


0 


0.1 1 0.205 


92 


25 


20.7 


7 


0 


0.1 


0.164 


2042 


68 


16.4 


22 


3 


0.9 


0.852 


2043 


63 


14.5 


9 


0 


0-9 


0.733 


2044 


84 


17.6 


23 


17 


0.9 


0.931 


2046 


72 


17.8 


8 


15 


0.9 


0.813 


2047 


61 


13.4 


146 


15 


0.9 


0.852 


2048 


44 


20.8 


21 


0 


0.9 


0.567 


2049 


60 


17.8 


29 


0 


0.9 


0.822 



An evaluation of the training data set {n=61) yielded a sensitivity of 96*^ 
and a specificity of 100%. 

Table II provides the data used to test the network trained with the 
5 data provided in Table I. 



TABLE n 

Testing Data for Ovarian Cancer Diagnostic System 



Patient 
Number 




LASA-Pe 


CA125 


DM/70K 


Histological 
Diagnosis | 


Neural Net 
Diagnostic 
Index 


103 


43 


14.1 


8 


15 




0.371 


104 


28 


15.9 


11 


0 


0.1 


0.211 


105 


46 


15.9 


15 


0 


0.1 


0.529 


107 


44 


23.6 


10 . 


26 


0.1 


0.416 


108 


32 


21.1 


9 


3 


0.1 


0.258 


109 


42 


14.1 


13 


12 


0.1 


0.407 


111 


35 


19.3 


9 


6 


0.1 


0.29 


113 


39 


18.7 


4 


7 


0.1 


0.305 


115 


31 


21.8 


10 


0 


0.1 


0.259 


2002 


69 


16.4 


29 


21 


0.9 


0.869 


2003 


53 


13 


11 


0 


0.9 


0.599 


2004 


58 


34.3 


80 


90 


0.9 


0.823 
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zUil 


o3 


14.3 


ID 


4 


A O 


A 017 

11.7 I / 


2012 


OU 


24.3 


2U 


3o 


A O 


U. / 33 


2013 


O 1 

81 


21 . / 


12 


lUO 


A O 


A fifi 


2014 


69 


33.4 


12 


39 


0.9 


A "700 


2016 


55 


17.3 


49 


5 


0.9 


0.821 


2017 


68 


19.6 


11 


0 


0.9 


0.813 


2018 


56 


17.8 


94 


3 


0.9 


0.832 


2021 


68 


12.2 


5 


0 


0.9 


0.761 


2023 


45 


15.9 


71 


1 


0.9 


0.713 


2025 


61 


17.9 


34 


5 


0.9 


0.846 


2028 


43 


17.6 


30 


0 


0.9 


0.617 


2029 


42 


16 


16 


0 


0.9 


0,47 


2030 


51 


21 J 


820 


57 


0.9 


A rr^ 
0.75 


47 


31 


21.7 


5 


0 


0. 1 


A TTJ 

U.223 


48 


39 


17.1 


3 


1 


O.l 


0.3U4 


50 


30 


17.2 


5 


2o 


A 1 

U. 1 


A liv 
U. lO 


54 


44 


21.1 


4 


0 


U.l 


A >l 1 1 

U.41 1 


55 


26 


t O f 

18.5 


1 


2 


A 1 

0. 1 


U. 133 


64 


23 


17.2 


2. 


u 


A 1 


A 1 1 


OD 




1 7.0 




1 

1. 


0 1 


0,483 


67 


27 


14.7 


3 


0 


0.1 


0.149 


69 


29 


17.2 


3 


12 


0.1 


0.158 


72 


26 


20.9 


2 


4 


0.1 


0.142 


74 


31 


19.8 


3 


20 


0.1 


0.173 


11 


25 


15 


10 


0 


0.1 


0.168 


78 


25 


12.4 


4 


0 


0.1 


0.132 


79 


35 


20,7 


1 


6 


0.1 


0.233 


81 


26 


17.2 


9 


0 


0.1 


0.178 


83 


45 


17.1 


2 


10 


0.1 


0.372 


86 


32 


15.6 


11 


0 


0.1 


0.262 
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87 


31 


21.5 


6 


11 


0.1 


0.21 


oo 


26 


17.2 


3 


0 


0.1 


0.145 






16.5 


2 


2 


0.1 


0.262 


91 


30 


16.9 


7 


0 


0.1 


0.211 


93 


34 


21.7 


8 


1 


0.1 


0.285 


95 


35 


20 


13 


0 


0.1 


0.339 


96 


34 


22.4 


7 


3 


0.1 


0.275 


97 


37 


15.9 


6 


6 


0.1 


0.284 


98 


38 


23.5 


38 


0 


0.1 


0.624 


2033 


72 


17.6 


3 


0 


0.9 


0.806 


2035 


38 


21.6 


19 


12 


0.9 


0.422 


2041 


59 


21.5 


1660 


9 


0.9 


0.857 


2045 


61 


23 


436 


46 


0.9 


0.844 


2050 


55 


28.4 


480 


139 


0.9 


0.794 



The testing data (n=59) yielded a sensitivity of 92% and specificity 
of 91%. 

The combined results of the training and testing data yielded a 
5 sensitivity of 94% and specificity of 96%. 

Fig. 20 provides the equations describing the weight matrices 
between the layers in the neural network 

Table ni below compares the sensitivity and specificity of a 
diagnosis obtained by the method of the present invention with the 
10 sensitivity and specificity of each individual biomaiker. As is evident from 

the results, the individual biomarkers may be combined by this method to 
yield a single diagnostic index having greater sensitivity than the individual 
tests for purposes of diagnosis. 
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TABLE m 

Comparision of Ovarian Cancer Diagnostic Methods 



BioMarker Test 


Sensitivity 


Specificity 


Carbohydrate Antigen 125 
(CA 125) 


74%* 


99%* 


Dianon Maiker 70K 
(DM/70K) 


68%* 


99%* 


Lqnd Associated Sialic Acid 
(LASA-P®) 


79%* 


96%* 


Neural Network Diagnostic 
Index for Ovarian Cancer 


94%** 


96%** 



* BioStatistics provided by Dianon Systems, Inc. 

BioSutisdcs based on initial pilot study performed by Horus 
Therapeutics, Inc. 

Example 7 



This example illustrates the construction and training of a neural 
network for diagnosis of osteoporosis. 

Fig. 9 provides the data used to develop, i.e. train and test, the 
neural network to diagnose osteoporosis. The biomarkers selected included 
age, calcium, phosphate, estradial (ETWO), progesterone, total alkaline 
phosphatase, total intestinal alkaline phosphatase, and % liver alkaline 
phosphatase. Fig. 9 further includes the diagnostic index obtained by the 
neural network. 

Fig. 10 provides the data use to test the netwodc trained with the data 
in Fig. 10, and the neural network diagnostic index obtained. 

The second embodiment of the present invention also relates to a 
method for diagnosing diseases in humans or animals using blood 
concentrations of specific, predetermined blood constituents. More 
specifically, the present invention relates to a method for detecting the 
occurrence of osteopenia, diagnosing osteoporosis and determining its 
severity and underlying cause. The present invention also facilitates the 
periodic monitoring of specific physiological functions which may indicate 
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the onset of osteopenia and conelates to bone mineral density measurements 
determined by various standard methods. 

In practicing one aspect of the present invention, the severity of 
disease in a set of humans or animals with varying severity of disease is 
measured by a standard method or methods. The measurement is then 
assigned a numerical value corresponding to a severity scale. The scale 
ranges from humans or animals with no disease, to humans or animals with 
severe disease. The scale is preferably a numerical scale. For example, one 
could assign a value which corresponds to normal or slight disease, another 
value which corresponds to moderate disease and a third value which 
corresponds to severe disease. 

The concentration of a predetermined set of blood constituents in the 
set of humans or animals with varying severity of disease is then 
determined. According to the present invention, it is preferable to measure 
the blood constituents in the same set of humans or animals in which the 
severity of disease was measured by the conventional method or methods. 

After determining the blood concentrations of the predetermined set 
of blood constituents in the set of humans or animals with varying severity 
of the disease, mathematical manipulations are performed in which the 
numerical value obtained using the conventional diagnosis is set to equal the 
result of a mathematical model utilizing the set of blood constituent 
measurements. For example, the relationship can be made using a multiple 
linear regreission analysis. It is to be understood that other mathematical 
models could be used to determine a correlation between the serum 
constituent concentrations and a standard method of determining severity of 
the disease. The concept of using mathematical models to detennine the 
relationship or correlation is considered to be part of the present invention. 
Standard statistical analyses that are well-known to tiiose of ordinary skill in 
the art are then performed to determined the confidence level of the 
correlation between the diagnosis by conventional means and the set of 
blood constituents. These statistical analyses include chi-square tests. 

An example of practicing one embodiment of the present invention is 
a method for diagnosing osteopenia in a human or animal. The method 
preferably utilizes six blood constituents. These constituents are calcium, 
phosphate, total alkaline phosphatase, an alkaline phosphatase isoenzyme, 
estradiol, and progesterone. The alkaline phosphatase isoenzymes preferred 
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for practicing the present invention include lymphocyte-derived alkaline 
phosphatase isoenzyme and bone, liver or intestinal allcaline phosphatase 
isoenzymes. The present invention includes calculating a bone density 
quotient using the aforementioned six blood constituents by entering the 
values for the tests into an algorithm that is calculated using a multiple linear 
regression analysis. Optionally, the age of the patient may also be 
incorporated into the equation. The bone density coefficient derived from 
the algorithm allows one to diagnose the osteopenic state of the patient, 
including the severity of the disease. 

In addition to diagnosing the osteopenic state of the human or 
animal, an indication of the underlying cause of the osteopenia can be 
determined using the present invention. For example, by practicing the 
present invention as described herein, one can determine whether the 
osteopenia in a human or animal is caused by post-menopausal lack of 
estrogen or is caused by some other condition, such as cancer. This allows 
the attending physician to be better able to prescribe the appropriate 
treatment for the osteopenia. 

Five of the serum tests that are used in the present invention are tests 
that are commonly performed by clinical laboratories. The test for t- 
lymphocyte derived alkaline phosphatase is experimental only; however, the 
test for blood, liver and intestinal alkaline phosphatase isoenzymes are also 
known. The type of test used to determine the six serum constituents is not 
critical to the present invention as long as the tests give accurate blood 
concentrations of the constituents being measured. 

Calcium & Phosphorus: 

Maintenance of calcium and phosphorus homeostasis involves the 
participation of 3 major organs: the small intestine, the kidney, and the 
skeleton, and is regulated by various hormones. Calcium enters the body 
through the diet and is absorbed into the circulation from the small intestine. 
Calcium absorption occurs by two processes, active transport and passive 
transport. Approximately 98% of the calcium and 85% of the phosphorus 
in the adult is present in the skeleton primarily as hydroxyapatite, which is a 
crystal lattice compound of calcium, phosphorus, and hydroxide. 
Remaining calcium is present in extra cellular fluid, some types of tissue, 
and skeletal muscle. Phosphorus is combined with lipids, proteins. 
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carbohydrates, and other organic substances. Of critical importance to 
calcium homeostasis is the fact that less than 1% of total skeletal reservoir of 
calcium is rapidly exchangeable with extracellular fluid. In addition to its 
obvious importance in skeletal mineralization, it is also vital for blood 
5 coagulation, neuromolecular condition, maintenance of normal tone and 

excitability of skeletal and cardiac muscle, and preservation of all membrane 
integrity and permeability, particularly in terms of sodium and potassium 
exchange. 

Free or ionized calcium accounts for 50% of total calcium. About 
10 5% of total calcium is complexed with a variety of anions, particularly 

phosphate and citrate. The remaining 45% of calcium is bound to plasma 
proteins. Both ionized calcium and calcium complexes are freely dialyzable. 
Addotic and alkalotic conditions adversely affect the ionized calcium level in 
the blood. In metabolic bone disease such as hyper- or 
15 hypoparathyroidism, Pagets of bone. Vitamin D deficiency, and renal 

osteodystrophy, calcium and/or phosphorus levels are significanUy altered. 

Serum Calcium Determination: 

In one method of practicing the present invention, the procedure for 

20 the determination of calcium is based on the interaction of the chromogenic 

agent o-cresolphthalein complexone (Sigma Diagnostics Calcium Agent, 
Sigma Chemical Co.. St Louis. MO ) which complexes with calcium cation 
in an alkaline medium to form a purple colored complex which has an 
absorbance maximum at 575 nm. The intensity of the color measured at 575 

25 nm is directly proportional to the calcium concentration in the given sample. 

The o-cresolphthalein complexone contains 8-hydroxyquinoline which 
prevents interference from magnesium ions. Robertson, W.G.. et al., 
"Calcium Measurements In Serum and Plasma - Total and Inonized." Crit. 
Rev. Clin, Lab. Sci.. Vol. 11, pg. 271. (1979); Sigma Diagnostics Product 

30 Insert: Calcium. Procedure No. 587. Reissued May 1989. The calcium 

level in serum is in the range of about 9.2 to 1 1 .0 mg/dl which is composed 
of 3 distinct fractions. 
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Specimen Collection: 

In one method of practicing the present invention, serum or 
heparinized plasma are suitable for the test. Anticoagulants other than 
heparin should not be used. 

5 

Procedure: 

Sigma Diagnostics Calcium Reagent is used with Abbott Spectrum 
High Performance Diagnostic System for the quantitation of calcium in 
serum. It is generally considered good laboratory practice to run a 
10 calibration, linearity, and a quality control prior to the determination of 

patient samples. Calcium levels up to 16 mgs % can be measured by this 
method. 

Serum Phosphorus Determination: 

15 In one method of practicing the present invention, the procedure for 

the determination of serum phosphorus is based on the interaction of 
inorganic phosphorus with ammonium molybdate in the presence of sulfuric 
acid. BarrMll S., etaL, "Osteoporosis: A Possible Autoimmune Etiology," 
AnrL of Clin. Lab. Set., Vol. 17, pg. 255, (1987). This reaction produces 

20 unreduced phosphomolybdate complex. The absorbance of this complex at 

340 nm is directly proportional to the inorganic phosphorus present in the 
given sample. (Sigma Diagnostics Product Insert: Phosphorus, Inorganic, 
Procedure No. 360-UV. Previous Revision December 1985. Reissued 
July (1988)). Most of the phosphorus in extracellular fluid is inorganic and 

25 in the range of 2.4 to 4.7 mg/dl. 



Specimen Collection: 

Serum or heparinized plasma are preferred for the test. 
Anticoagulants other than heparin should not be used. 
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Procedure: 

Sigma Chemical Company's phosphorus inorganic reagent can be 
used with the Abbott Spectrum High Performance Diagnostic System for 
qualification of the inorganic phosphorus in serum. It is customary to 
5 perform calibration, linearity, and quaUty control studies prior to the 

determination of test samples. Phosphorus levels up to 12 mgs % can be 
measured by this procedure. 

Total Alkaline Phosphatase: 

IQ The demonstration that bone is rich in alkaline phosphatase (ALP) 

and that normal serum contains the same or a similar enzyme led to the study 
of serum ALP levels in patients with disease of bone. Sec. Courpron P.. 
"Bone Tissue Mechanisms Underlying Osteoporosis." Orthop, Clin. North 
Am., Vol. 12. page 513, (1981). This is especially true in osteitis 

15 deformans, hyperparathyroidism and bone neoplasm. ALP is also elevated 

in hepatic disease; however, this can be distinguished by other conoborative 
laboratory procedures and clinical feamres. In some situations, as in 
osteoporosis. ALP per se may not be above the reference range; however, 
the isoenzyme of ALP is increased. It is weU understood that total serum 

20 ALP in normal subjects consists of isoenzymes contributed from liver, 

bone, renal, pulmonary, placental and intestinal sources, among others. 
The isoenzymes and isoforms of various tissue origins can be further 
separated and visualized by an Isoelectric Focusing Technique. Estimation 
of total serum ALP by isoelectric focusing techniques offers great potential 
25 in the investigation of metabolic bone diseases. 

Serum Alkaline Phosphatase Determinanon: 

In one method of practicing the present invention, serum ALP 

activity can be measured using various phosphate esters as substrates. 
30 Sigma Chemical Company's alkaline phosphatase reagent measures serum 

ALP activity by a kinetic method. The reagent for the test contains p- 

nitrophenyl phosphate, carbonate buffer, magnesium ions and mannitol. 

Mannitol present in the reagent acts as a phosphate acceptor during the 

enzyme reaction. McComb iLb., et al. AlKRllTlfi PhOSffhaiase, Plenum 
35 Press. New York, (1979); Gundbery, M.. AlKRlinn Phosphatase and 

OQtencalcin. Primer on M ^fahnlic Bone DisftaSRS and DiSOrdfirS Of Minera l 
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Metabolism 1st Edition, Published by Am. Soc. For Bone and Mineral Res. 
Editor-Murry J. Favus, pgs. 74-76, (1990); Sigma Diagnostics Product 
Insert: Alkaline Phosphatase (ALP), Procedure No. 245, Reissued March 
1989; Kaplan T^.,et aL, Clinical Chemistry. Sl Louis, e.v. Mosby 
5 Company, (1987). 

Specimen Collection: 

Preferably, senun or heparinized plasma is used for the total alkaline 
phosphatase and alkaline phosphatase isoenzyme determination tests. 
ID EDTA, oxalate, citrate and fluoride are inhibitors of ALP and are not 

suitable anticoagulants. 

Procedure: 

Alkaline phosphatase (ALP) diagnostic reagent manufactured by 
15 Sigma Chemical Company is used with the Abbott Spectrum High 

Performance Diagnostic System for the quantitation of inorganic 
phosphorus in the serum. Serum ALP hydrolyzes p-nitrophenyl phosphate 
to p-nitrophenol and inorganic phosphate. The hydrolysis occurs at alkaline 
pH and the absorbance is directly proportional to the ALP activity of the 
20 serum sample. As is the usual practice, calibration, linearity and quality 

control assays should be performed prior to the determination of test 
samples. 

ALP levels up to 1200 U/L can be measured by this procedure. 
Normal ran^ are: 
25 Infants 50-165 (U/L) 

Adult 20-70 (U/L) 

Child 20-150 (U/L) 

60 years 30-75 (U/L) 

30 Alkaline Phosphatase Isoenzyme: 

Alkaline phosphatase isoenzyme levels have been primarily used to 
aid in the differential diagnosis of liver and bone disorders and also to 
indicate placental growth patterns. Isoelectric focusing of alkaline 
phosphatase isoenzyme on a precast gel gives 12 discernible bands of 

35 different organ and tissue origin. Of these 12 Bands, Band 10, appearing at 

4.73 isoelectric point, is of placental origin. It is released by placental 
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syncytiouophoblast and/or macrophage* and is T-lymphocytc mediated. In 
the alternative, samples for ALP isoenzyme can be analyzed by 
electrophoresis methods which identify mainly liver, bone, and intestinal 
isoenzymes. The Ciba Coming Alkaline Phosphatase Isoenzyme (ALP) 
5 System may be used for the qualitative and quantitative determination of 

ALP isoenzymes in human serum by electrophoresis. 

The liver isoenzyme is the most frequendy encountered alkaline 
phosphatase isoenzyme in pathologically increased serum alkaline 
phosphatase activity. Increased Uver isoenzyme is encountered in a variety 
10 of Uver and hepatobiliary diseases. These include hepatic cirrhosis, primary 

biliary cirrhosis, congestive cirrhosis, and hepatic carcinoma. Elevation of 
the liver isoenzyme is also a sensitive indicator of cholestasis and liver 
infiltration. 

Another of the alkaline phosphatase isoenzymes that can be used in 
15 practicing the present invention is lymphocyte-derived alkaline phosphatase 

The preferred method of measuring this isoenzyme is by isoelectric focusing 
electrophoresis. A commercially available isoelectric focusing apparatus thai 
is capable of measuring alkaline phosphatase isoenzymes is made by 
IsoLab, Inc. (Resolve® -ALP, IsoLab, Inc., Akron, Ohio). On the 
20 Resolve®-ALP isoelecuic focusing apparatus, the alkaline phosphatase 

isoenzyme that is used in the present invention resolves as band 10 in the 
elecuophoreUc pattern. In the preferred bone density algorithm defined 
hereinbelow, the value for the band 10 alkaline phosphatase isoenzyme was 
assigned 0 if the band were missing or very weak, and 1 if the band were 
25 present. The electrophoresis gels can also be scanned with a densitometer 

and more quantitative values can be assigned to the alkaline phosphatase 
isoenzyme concenttation. As used herein, the band 10 alkaline phosphatase 
isoenzyme shall be designated lymphocyte-derived alkaline phosphause 
isoenzyme, or I-Alkp. 

30 It is to be understood that other methods of measuring I- Alkp could 

be used. These methods include, but are not limited to, enzyme-linked 
immunoassay techniques (ELISA), radioimmunoassay techniques, affinity 
columns, and isoelectric focusing columns. It is important to note that 
measurement of band 10 on the Resolve®-ALP apparatus does not correlate 

35 at aU with osteopenia in general There are many other abnormal conditions 
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which result in a higher than normal band 10 lymphocyte-derived alkaline 
phosphatase isoenzyme (see Bamhill, et al. and Griffith, et aL^ supra.). 

A number of pathological conditions can lead to an elevated level of 
bone alkaline phosphatase isoenzyme. The highest levels of bone 
isoenzymes are usually found in Pagets disease. Increased levels are also 
encoimtered in rickets, bone cancer, osteomacacia and celiac sprue. Renal 
disorders can also result in increased levels. These include renal failure, 
primary hyperparathyroidism, secondary hyperparathyroidism induced by 
long term hemodialysis and malabsorption. 

Increased levels of intestinal alkaUne phosphatase are encountered in 
a variety of diseases of the digestive tract These include intestinal infection 
and ulcerative lesions of the stomach, duodenum, small intestine and colon. 

Procedure: 

In the preferred method of practicing the present invention, the 
alkaline phosphatase isoenzymes are separated by electrophoresis in a 
buffered agarose system. After electrophoresis the isoenzymes are deu:cicd 
by incubating the gel with a fluorescent compound such as 4- 
methylimibelliferyl phosphate. The fluorescence formed during the 
interaction is quantitated using a Ciba-coming 710 densitometer at 385 nm. 
Total alkaline phosphatase activity of the patient is required to interpret the 
liver, bone and intestinal ALP isoenzymes. 

Estrogen and Progesterone: 

The estrogens are steroids that have a ring containing three 
unsamrated double bonds. The ovary, as well as the testes and adrenal 
gland, has the capacity to synthesize estrogens from androgens, 
androstenedione and testosterone. During the follicular phase of the 
menstrual cycle, ovarian secretion represents only one third of total estrogen 
production. In contrast to estradiol, which is secreted almost entirely by the 
ovary, most estrone is derived from peripheral conversion of 
androstenedione and from estradiol metabolism. During menopause, 
estradiol concentrations steadily decrease to approximately 15 percent of 
premenopausal levels. In healthy postmenopausal women, the ovaries do 
not secrete significant quantities of estrogens, and virtually all estrogen 
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produced is from peripheral conversion of androstenedione made by the 
adrenal. Estradiol is present in the serum as follows: 



Follicular Phase 




Early 


30-100 ng/L 


Late 


100-400 np/L 


Luteal Phase 


50-150 ne/L 


Postmenapausal 


<20ne/L 



10 



15 



In menstruating females, progesterone is secreted mainly by the 
corpus luteum of the ovary. It is partially responsible for cyclic changes in 
the endometriiun that are necessary for attachment and growth of an 
embryo. Progesterone levels are low prior to the mid-cycle gonadotropin 
surge. Shortly after the gonadotropin surge, they begin to rise rapidly, 
reaching peak levels during the middle of the luteal phase. Thereafter, a 
progressive fall occurs with barely detectable progesterone levels prior to 
menses. Although progesterone in large amounts produces a negative 
feedback on gonadotropin secretion, it is not the major component in the 
negative feedback system of ovarian steroids. Function of the corpus 
luteum can be assessed by measuring serum progesterone concentration. 
Progesterone is present in the serum as follows: 



Follicular Phase 


0.1 - 1.5 ng/L 


Luteal Phase 


2.5 - 28.1 ng/L 


Mid-Luteal Phase 


5.7 - 28.1 ng 


Over 60 Years 


0.0 - 0.2 ng/L 



20 



Procedure: 

In a preferred method of practicing the present invention, progesterone 
measurements can be obtained by a radioimmuno-assay method using 
antibody coated tubes (Diagnostic Products, Los Angles, California). 
Estradiol measurements can be performed by a Microparticle Enzyme 
Immunoassay (IMX). available from Abbott Diagnostics, Abbott Park, IL. 



25 
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Bone Mineral Density (BMD) Measurements 

The conventional methods of diagnosing osteopenia which may be 
used when practicing the present invention are outlined below: 

Dual-Photon absorptometry (DP A) is widely used to assess bone 
mineral content and bone mineral density. Johnston Conrad C, et al^ 
**Clinical Use of Bone Densitometry'*, New Eng. J. of Med., Vol. 324, pgs. 
1 lOS-1 109, (1991). DPA uses transmission scanning with photons from a 
radioisotope source, such as 153 Gd, that emits two energy peaks, thus 
allowing bone density to be measured independent of soft tissue. DPA 
measurements are performed on lumbar spine (L1-L4) and femur (femoral 
neck. Ward's Triangle and trochanter) and the average determined 
separately. The overall average of both hip and spine can also be 
determined. BMD can be measured by DPA using - Gadolinium as the 
source (Lunar DP3, by Lunar Radiation Corporation, Madison, Wisconsin). 

The present method employs transmission scanning using 44 and 
100 KeV photon energies from a one Ci^^^Gd source to allow computation 
of the mineral content of bone independent of soft tissue thickness. Bone 
mineral density, expressed in g/cm2, is derived by dividing bone mineral 
content (BMC) by the projected area of the scanned bone. Peppier W,W., 
et al, "Total Body Bone Mineral and Lean Body Mass by Dual-Photon 
Absorptometry: 1. Theory and Measurement Procedure," CALCIP, Issue 
Int., Vol. 33, pg. 353, (1981); Shipp C, etoL, "Precision of Dual-Photon 
Absorptometry." CALCIP, Issue Int., Vol. 42, pgs. 287-292, (1988). 
During spine scans, the detector moves in a rectilinear pattern at a rate of 
5mm/sec and with scan lines 4.5 mm apart The BMC and BMD are 
calculated in lumbar vertebrae 1 through 4 (including intervertebral discs) 
with a software version supplied by Limar Radiation Corporation of 
Madison, Wisconsin. 

During femur measurements, the scanner moves at a rate of 2.5 
mm/sec and a step distance of 2.5 mm. The BMC and BMD of the neck. 
Ward's Triangle and trochanteric regions of the proximal femur are 
calculated using a femur software version supplied by Lunar Radiation 
Corporation. The femoral neck region of interest (ROI) is that band about 
1.5 cm wide across the neck of the bone perpendicular to the neutral axis 
with the lowest density. Ward's Triangle is defined as a square ROI (about 
1.5 X 1.5 cm) with the lowest density within the proximal femur region. 
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Ward's Triangle is predominantly trabecular bone and contains the least 
amount of bone mineral within the neck region. Carter D.R., et al., 
"Relationship Between Loading History and Femoral Cancellous Bone 
Architecture/ J. Biomechanics. Vol. 22, pgs. 231-244, (1989). In the 
proximal femur, the ROI is usually the area 1.5 cm wide across the entire 
femoral neck. Additional regions are defined by the software in the lower 
density Ward's Triangle region, and in the region of the greater trochanter. 
Bone loss in the proximal femur begins in the Ward's Triangle region and 
proceeds outward from there. {Brown £>.. et al,, "Mechanical Property 
Distributions in the Cancellous Bone of the Human Proximal Femur," Act. 
Orthop. Scand.. Vol. 51, pgs. 429^37, (1990). This makes the region an 
early indicator of bone loss, but the higher variance in measuring it, 
compared to the neck region, makes the latter zone, a better discriminator. 
However, Ward's Triangle is the lowest density area at the point where the 
neck and greater trochanter meet, a primary hip fracture site. This 
operational ROI may not correspond exactly to the anatomic Ward's 
Triangle region but does provide a repeatable measurement. The width of 
the neck ROI and the size of the Ward's Triangle ROI are actually 
proportional to the measured size of the femoral neck. It has been shown 
that the density of the Ward's Triangle area is substantially reduced in hip 
fracture patients compared with age matched controls. {Vose G. P., et oL, 
"Femoral Neck Fracturing its Relationship to Radiographic Bone Density," 
J. Gerontol, Vol. 20, pgs. 300-305, (1965). DPA measurements have a 
precision of 1-3% and the scan can be completed in about 20 minutes. Bone 
density of vertebrae correlates well with risk for vertebral fracture, and bone 
density of areas in the proximal femur correlates well with risk for hip 
fracmre. DPA has a low radiation dose (<10 mrem to skin, and 2 mrem to 
marrow). DPA has approximate sensitivity of spine and hip of 50% and 
53%. respectively, at 95% specificity (i.e. % of fracture cases below 5th 
pcicentae). Moreover. DPA is not subject to systematic errors introduced 
by variable osteoid and variable marrow. 

Bone density can also be measured by radiographic absorptomeury. 
Radiographic absorptometry is a method of measuring bone density which 
is well known to those of ordinary skill in the an. Other methods of 
measuring bone density include quantitative computed tomography and 
direct measurement of bone density. 
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Direct measurement of bone density can be obtained by measuring 
the bone density in cadavers. Thus, another method of determining the 
correlation between the blood concentrations of certain predetennined blood 
constituents and bone density is to measure the bone density of a set of 
cadavers from which blood concentrations of predetermined blood 
constituents are known. The multiple linear regression analysis can then be 
performed and a bone density coefficient relationship can easily be 
determined. 

In correlating the bone density measurements to concentration of 
blood constituents, blood concentrations of calcium, phosphate, total 
alkaline phosphatase, an alkaline phosphatase isoenzyme estradiol, and 
progesterone are measured. Liver, bone and intestinal isoenzymes can be 
used. Band 10 alkaline phosphatase isoenzyme. A mathematical 
relationship between the concentrations of blood constituents and bone 
density, as measured by radiographic absorptometry or other standard 
method of measuring bone density, is determined by performing a multiple 
regression analysis with the following model: 

Bone density coefficient =*o + ^^iCa + b2P + fc3E2 + 64Pg +fe5Alkp + 
fe6l-Alkp + &7(CaCa) + bi(PP) + fe9(E2E2) + felO(Pg + Pg) + 
fell(AlkpAlkp) + frl2(CaP) + fci3(CaE2) + £7i4(CaPg) +fei5(CaAlkp) + 
fcl6(CaI-Alkp) + fcl7(PE2) + fel8(PPg) + frl9(PAlkp) + fc20(PI-Alkp) + 
&21(E2Pg) + £?22(E2Alkp) + fc23(E2I-Alkp) + i^24(PgAlkp) + 625(Pg^ 
Alkp) -h 2;26(AlkpI-Alkp) 
where: 

Ca = Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 

I- Alkp = Serum concentration of alkaline phosphatase isoenzyme. 

Calculation of the mathematical model utilized the Systat® statistical 
package (Systat: Inc., Evanston, IL). The multiple linear regression 
analysis is an iterative process which calculates the correct coefficients so 
that the result of the algorithm using the blood concentrations correlates to a 
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high degree with the result of the bone density measurement by radiographic 
absorptometry. 

For diagnosing osteopenia, the general form of the preferred 
algorithm that is used in the present invention is as follows: 

Bone density coefficient = fco + MCa + fr2P + i^E2 + i>4Pg -i-^^sAlkp + 
f>6l-Alkp + d7(CaCa) + t8(PP) + fe9(E2E2) + ilO(Pg + Pg) + 
fcll(AlkpAlkp) + fel2(CaP) + fei3(CaE2) + &l4(CaPg) +£>i5(CaAlkp) + 
frl6(CaI-Alkp) + tl7(PE2) + fcl8(PPg) + fcl9(PAlkp) + &20(PI-Alkp) + 
fe2l(E2Pg) + i^2(E2Alkp) + t23(E2I-Alkp) + b24(PgAlkp) + fe25(Pgl- 
Alkp) + &26{AlkpI-Alkp) 
where: 

bo is a constant 

Ca - Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Senmi concentration of total alkaline phosphatase 

I-Allqp = Serum concentration of alkaline phosphatase isoenzyme. 

It is to be understood that concentration of I- Alkp can be replaced 
with serum concentration of liver, bone or intestinal derived alkaline 
phosphatase isoenzymes. 

In practicing the present invention, the values obtained from the 
patient for the six blood constituent concentrations are inserted into the 
algorithm as indicated. The mathematical manipulation is performed. The 
resulting number is called a bone density coefficiait and is then placed in the 
severity scale. This results in a bone density probabiUty quotient which 
correlates to a high degree to bone density as measured by radiographic 
methods. 

As seen in the algorithm that is considered part of the present 
invention, there are a number of coefficients which are a part of the 
algorithm. It is to be understood that these coefficients can change if the 
bone density coefficient is correlated to a different method of determining 
bone density. Thus, if bone density is measured using dual photon 
absorptometry, and one wanted to correlate the bone density coefficient to 
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the results of the dual photon absorptometry, the overall relationship of the 
tests in the algorithm would be the same or similar as disclosed herein but 
the coefficients could be diffenmt 

In its preferred embodiment, the bone density probability quotient is 
assigned to one of the following diagnostic categories: 

Normal to mild osteopenia 
Modoate osteopenia 
Severe osteopenia 

In addition, the menopausal status may be determined based on the 
results of the estradiol and progesterone levels into one of the following 
diagnostic categories: 

Probable Pre-menopausal 
Probable Peri-mmopausal 
Probable Post-menopausal 

Thus, according to the present invention, using blood concentrations 
of certain blood constituents, one can not only diagnose bone density but 
can also obtain an indication of the underlying cause of the osteopenia. 

This invention is further illustrated by the following examples, 
which are not to be construed in any way as imposing limitations upon the 
scope thereof. To the contrary, it is to be clearly imderstood that resoa may 
be had to various other embodiments, modifications, and equivalents thereof 
which, after reading die description herein, may suggest themselves to those 
skilled in the ait without departing from the spirit of the present invention 
and/or the scope of the appended claims. 

Example 8 

In the following example, the bone density coefficient algorithm was 
correlated with bone density as measured by radiogr^hic absorptometry. 
The coefficients were calculated by conelating the blood concentration of the 
indicated blood constituents in 27 patients to the results of bone density in 
those patients as measured by radiographic absorptometry. The bone 
density algorithm which results is as follows: 
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0 + (-8.955)(Ca) + (79.370)(P) + (34.076)(E2) + (-9.216)(Pg) + 
(.0.600)(Alkp) + (-57.855)(I-Alkp) + (0.926)(CaCa) + (-2.735)(PP) + 
(.0.272)(E2E2) -h (-0.064)(PgPg) + {0.0{)4)(AlkpAlkp) + (-4.029)(CaP) + 
(-3.156)(CaE2) + (1.172)(CaPg) + (-0.031)(CaAlkp) + (8.238)(CaI-Alkp) 
+ (1.906)(PE2) + (-0.138)(PPg) + (-0.018)(PAlkp) + (-6.710)(PI-Alkp) + 
(0.193)(E2Pg) + (-0.048)(E2Alkp) + (.1.825)(E2I.Alkp) 

wherein: 

Ca = Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of serum lymphycyte-derived alkaline 

phosphatase isoenzyme. 

The bone density coefficient that is obtained from the bone density 
algorithm is assigned to one of three categories: 

Group 1 normal to mild osteopenia 

Group II moderate osteopenia 

Group m severe osteopenia 

Group in ^ approximately 78.5 < Group n ^ approximately 100 < Group I 
It has been determined that the bone density coefficioit calculated according 
to the present invention places patients in the same groups as radiographic 
absorptometry measurements for bone density. Thus, the present invention 
provides a safe, economical and accurate method for diagnosing osteopenia. 

Bone density quotients derived practicing the present invention using 
the six biochwnical serum constituents were found to strongly correlate with 
bone density measurements from radiographic absorptometry. It is to be 
understood that the preferred embodiment utilizes all six of the biochemical 
serum constituents listed hereinabove. However, if any five of the serum 
constiments are utilized and then included in the bone density algorithm, a 
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bone density quotient is obtained which does not conelate as weU as when 
six tests are used. 

It is contemplated as part of the present invention that any five of the 
biochemical senim constituents can be used to determine the osteopenic sute 
of a patient For example, if senim calcium is removed from the bone 
density algorithm, then the algorithm predictive value for moderate or severe 
osteopenia is approximately 66% altiiough tiie predictive value for normal to 
mUd osteopenia is 100%. If progesterone is removed from die bone density 
algorithm, then the algorithm predictive value for moderate or severe 
osteopenia of approximately 75% and a predictive value for normal to mild 
osteopenia of approximately 94%. If alkaline phosphatase is removed from 
die bone density algoritiun. tiien tiie algoritiim predictive value for moderate 
or severe osteopenia of approximately 83% and a predictive value for 
normal to mUd osteopenia of approximately 85%. If estiadiol is removed 
from the bone density algoritiun, tiien the algoritiim predictive value for 
moderate or severe osteopenia of approximately 67% and a predictive value 
for normal to mild osteopenia of approximately 85%. However, it has been 
determined tiiat when aU six tests are included in the bone density algoritiun. 
tiie predictive value for moderate and severe osteopenia is 100% and tiie 
predictive value for normal to mild osteopenia is 100%. AU of tiie predictive 
values may vary sUghtiy witii dififerent patient populations. 



25 woman: 



Example 9 

The foUowing serum constituents are measured in a 37 year old 



Calcium 9.5mg/dL 

Phosphate 5.0mg/dL 

Estradiol 2pg/mL 

Progesterone 0.2ng/mL 

Alkaline Phosphatase go U/L 

Alkaline phosphatase isoenzyme o (negative) 
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The above blood chemisuy values were used to calculate a bone 
density coefficient using die following algoritiun: 
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0 + (-8.955)(Ca) + (79.370)(P) + (34.076)(E2) + (-9.216)(Pg) + 
(-a600)(Alkp) + (-57.855)(I.Alkp) + (0.926)(CaCa) + (-2.735)(PP) -i^ 
(-0.272)(E2E2) + (-0.064)(PgPg) + (0.004)(AlkpAlkp) + H.029)(CaP) + 
(-3.156)(CaE2) + (1.172){CaPg) + (-a031)(CaAlkp) + (8.238)(CaI-Alkp) 
5 + (1.906)(PE2) + (.0.138)(PPg) + (-0.018)(PAlkp) + (-6.710)(PI-AUcp) + 

(0.193)(E2Pg) + (.a048)(E2Alkp) + M.825)(E2I.Alkp) 
where: 

Ca = Serum concentration of calcium 
P = Serum concentration of phosphorus 
10 E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of serum lymphycyte-derived alkaline 

phosphatase isoenzyme. 



A bone density coefficient of 101.2 is calculated for this paUent. 
This coefficient falls within Group I in the following severity scale: 

Group m ^ 78.5 < Group II ^ 100 < Group I 

Group I normal to mild osteopenia 

Group II moderate osteopenia 

Group ni severe osteopenia 

This patient has normal to mild osteopenia. When the patient's bone 
density is measured by radiographic absorptometry, she is found to have 
normal bone density. 
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Example 10 

The following serum constiments are measured in a 47 year old 
woman: 
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Calcium , 9.5mg/dL 

Phosphate 5.0mg/dL 

Estradiol 2pg/niL 

Progesterone o.2 ng/ml 

Alkaline Phosphatase go U/L 

Alkaline phosphatase isoenzyme l (positive) 



The values are inserted into the bone density algorithm recited in 
Example I and a bone density coefficient of 84.4 is calculated. This falls in 
Group n, indicating a moderate osteopenia. When the patient's bone 
density was measured by radiographic absorptometry, she was found to be 
moderately osteopenic. 

Example 11 

The following serum constituents are measured in a 47 year old 
woman who has been diagnosed with breast cancer: 



Calcium 9.0mg/dL 

Phosphate 4.0 mg/dL 

Estradiol Opg/mL 

Progesterone 0 ng/mL 

Alkaline Phosphatase 60 U/L 

Alkaline phosphatase isoenzyme 1 (posidve) 



The values are inserted into the bone density algorithm recited in 
Example I and a bone density coefficient of 69.9 is calculated. This falls in 
Group m indicating a severe osteopenia insured by another method. 

Example 12 

The foUowing example shows the correlation between the method of 
the present invention and Dual Photon Absorptometiy (DPA) as a tool for 
diagnosing and determining the severity of osteoporosis in a particular 
individual. 200 female subjects representing a cross section of ages and 
menopausal status were evaluated. These specific age groups were 
represented as follows: 
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As NumhernfP^it^fTir^ 

25-35 26 

36-45 25 

46-55 33 

56-65 30 

66-75 34 

76-85 25 

85 and over 26 



Each subject completed an Osteoporosis Data Questionnaire Form, 
which included detailed information on personal history, family history, 
gynecologic history, medical history, surgical history, and drug history, 
with special references to hormone treatment. 

Three tubes of blood were collected from each subject via 
venipuncture. One tube was sent to Bamhill-MetPath Laboratories in 
Savannah, Georgia for analysis of serum calcium, phosphorous, and total 
alkaline phosphatase. The second tube was sent to Horus Therapeutics, 
Inc. in Savannah, Georgia for determination of lymphocyte alkaline 
phosphatase utilizing the Isoelectric Focusing Method. The third tube was 
sent to MetPath Laboratories in Teterboro, New Jersey for evaluation of 
estradiol and progesterone by immunoassay. Data from all three sources 
were then subjected to statistical analysis. At the time of the analysis of 
data, each subject's chart was critically reviewed for relevant clinical 
information. 

The Bone Mineral Density (HMD) of each subject was also 
measured by Dual-Photon Absorptometry using 153 Gadolinium as a source 
(using Lunar DP3, Lunar Radiation Corporation, Madison, Wisconsin). 
DPA measurements were performed on the Lumbar Spine (L1-L4) and the 
hip, specifically the femoral neck. Ward's Triangle, and the trochanteric 
regions. HMD was measured in gm/cm^ for each region noted, and the 
fracture risk was determined. 

The BMD measurements were then age matched, and adjusted for 
sex, age, ethnic group, and weight. From the individual values obtained, 
average BMD for hip and spine was calculated. The overall average BMD 
for both hip and spine was also computed. 
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A mathematical algorithm correlating the results using the present 
invention with DPA was determined with the resultant data is shown in 
Figure 11. 

Figure 12 illustrates the discrepancies between the DPA 
5 measurements and the measurements obtained using the method of the 

present invention for individual subjects. The amount of the discrepancy 
(the residual) for any patient is the vertical distance between the points for 
that patient on the DPA and the present invention lines. 

Figure 13 illustrates the appropriateness of using the present 
10 invention as a prediction of DPA measurements. Note the absence of any 

tendency for points to be consistently above or consistently below the 
diagonal straight line* as greater DPA measurements are considered. 
Additionally, there is an absence of any tendency for points along any 
vertical line to be more widespread than are points along any other vertical 
15 line. Finally, note that the generally elongated shape foraied by the points is 

consistent with a strong linear correlation between DPA and the present 
invention. 

Figure 14 demonstrates the normality of values predicted by the 
present invention (the horizontal scale). The vertical scale shows the of the 

20 normalized scale (z) scores values determined according to the present 

invention as would be expected from the assumption of the present 
invention values fom[iing a normally distributed collection of numbers. 

The closer the dots are to the straight line, the greater is the 
consistency between what is observed about the present invention values 

2S and what is expected about them from the assumption of their forming a 

nomally distributed collection of numbers. 

Figure 15 demonstrates the normality of residuals (the difference 
between the observed DPA measurement and the value predicated by present 
invention), as evidenced by the arrangements of points along the diagonal 

30 straight line. The horizontal scale shows values of residuals The vertical 

scale shows the residuals' normalized scale (z) scores as would be expected 
from the assumption that the residuals form a normally distributed collection 
of ntambers. 

The closer the dots are to the straight line, the greater is the 
3S consistency between what is observed about the residuals and what is 
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expected about them from the assumption of their forming a normally 
distributed collection of numbers. 

Figure 16 demonstrates that the variance of residuals stays the same 
regardless of the value calculated according to the present invention involved 
in the residuals' computation. This homogeneity of variance appears in the 
uniform density of the dots in a rectangle formed by residuals between -0.17 
and 40.17 and by the present invention between 40.4 and +1.0. Moreover, 
there is no discemable tendency for points along any vertical line to be more 
widespread than are points along any other vertical line. 

Figure 17 demonstrates the independence of residuals and the values 
predicted by the present invention. The horizontal scales shows the values 
predicted by the present invention. The vertical scale shows the residual's 
studentized scale (t) scores. 

When there is no relationship between values predicted by the 
present invention and the residuals involving those predicted values (that is, 
when there is independence between residuals and predicted present 
invention values) the following criteria should be met: 

1 . There should be no discemable pattern (such as a straight or 
curved line) to the dots. 

2. About half the dots should appear in the top half of the 
graph. 

Figure 18 demonstrates that virtually all observed DPA 
measurements are satisfactorily predicted by the method of the present 
invention. The horizontal scale shows values predicted by the method of the 
present invention. Ilie vertical scale shows Cook's distances between 
observed DPA measurements and the predictions of those measurements by 
the method of the present invention Thus, the closer the dots are to a height 
of Cook = 0, the more accurate is the method of the present invention. 
Because Cook's distances tend to form a collection of numbers fitting an F- 
statistical distribution, the two outline points at heights between 10 and 20 
are not thought to be problematic in view of the large size of the sample. 
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Example 13 

The following example describes the training of a neural network to 
prognose prostate cancer. 

As shown in Fig. 19, a total of 52 samples were divided into 2 

5 groups* a training set and a generalization testing set. The training set 

contained 40 samples (28 stable and 12 progressing) and the generalization 
testing set contained 12 samples (9 stable and 3 progressing). 

The initial network architecture was selected based on the level of 
complexity of the classification task. A multi-layer feedforward network 

10 was used. Selection of the initial architecture involved the selection of the 

number of hidden layers and the nimiber of neurons in each hidden layer. 
Several trial iterations were performed to determine and adequate 
configuration that showed good results on both the training sample set and 
the generalization test sample set The present network had one hidden 

15 layer, having nine neurons, and two output neurons. 

Initially, connection weights among the neurons were randomly set 
The neural network had Ave input neurons, corresponding to five input 
variables significant for prostate cancer: TPS. PSA. PAP. CEA. and 
testosterone. The training data is shown in Fig. 19. During training, the 

20 five input variables for each patient were first linearly scaled into the 

continuos range between 0.0 and 1.0. The resultant five numbers were then 
presented as an input vector to the input neurons of the artificial neural 
network. 

For each of the input vectors, the netwoik generated an output based 
25 on the connection weights among the netwoik neurons. The output can be a 

single value or a vector of numbers, depending on the number of output 
neurons used. The network used had two output neurons. The outputs of 
the two neurons were processed by the following mathematical equation to 
yield a single diagnostic index: 



30 



Index = fANN2-ANNn + 0.5 
2 



35 



Each neuron in the network participated in the output calculation by passing 
the sum of all inputs to the neuron through a non-linear s-shaped function 
(often a ligistic function) and sending the result to each and everyone of the 
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neurons in the following adjacent layer. The generated output or each 
output neuron was compared to the desired "target** output. A value of 0.1 
corresponded to a diagnosis of stable and an output of 0.9 corresponded to a 
diagnosis of progressing. The difference was used to calculate an error term 
to guide the training algorithm, i.e., the back propagation algorithm, in the 
adjustment of network connection weights in an attempt to reduce the 
differences between network outputs and target values over the training 
sample set 

After training, the neural network correctly classified 100% of the 
samples. 

When presented with the generalization test results, the trained 
neural network correcdy identified 100% of the stable samples and 66% of 
the samples where the disease was progressing. 

It should be understood, of course, that the foregoing relates only to 
preferred embodiments of the present invention and that numerous 
modifications or alterations may be made therein without departing from the 
spirit and the scope of the invention as set forth in the appended claims. 
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Claims 

What is claimed is: 

1 . A method for diagnosing a disease in a human or animal 
comprising the steps of: 

measuring the concentrations of a predetermined set of biomaikers 
known to be associated with the disease from a biological fluid from 
the human or animal; 

scaling the digitized values of the biomarkers concentrations; and 

introducing the scaled values to a trained neural network means, 
whereby the output values from the neural network means tend 
toward the upper value when the human or animal has the disease 
and the output values tend toward the lower value when the human 
or animal does not have the disease. 

2 . An apparatus for diagnosing a disease in a human or animal 
comprising: 

means for digitizing the concentrations of a predetermined set of 
biomarkers known to be associated with the disease from a 
biological fluid from the human or animal; 

a means for scaling the digitized values; and 

a trained neural network means coupled to the digitizing and scaling 
means for generating network output values between an upper and 
lower value» whereby the ou^ut values tend toward the upper value 
when the human or animal has the disease and the output values trad 
toward the lower value when the human or animal does not have the 
disease. 
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3. A method for diagnosing a disease in a human or animal 
comprising the steps of: 

measuring the concentrations of a predetermined set of biomarkers 
known to be associated with the disease from a biological fluid from 
the human or animal; 

scaling the digitized values of the biomarkers concentrations; 

introducing the scaled values to a first trained neural network means, 
whereby the output values from the first neural network means tend 
toward the upper value when the human or animal has the disease 
and the output values tend toward the lower value when the human 
or animal does not have the disease; and 

introducing the output value from the first neural network and a 
second set of predetermined biomarkers to a second trained neural 
network means, whereby the output values from the second neural 
network means tend toward the upper value when the human or 
animal has the disease and the output values from the second neural 
network means tend toward the lower value when the human or 
animal does not have the disease. 

4. The method of Claim 3 wherein the second set of 
predetermined biomarkers includes patient biographical information. 

5. The method of Claim 3 wherein the second set of 
predetermined biomarkers are the same as the first set of predetermined 
biomarkers. 
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PATIENT 




CAL 


PH05 


ETWO 


PROG 


TALP 


TINTES 




SCORE 


143 


66 


9.4 


3.3 


14 


0«1 


71 . 


\J 








52 


9 


?<5 


5? 


0.3 


64 








145 


82 


8.5 


3.1 


■ c 


0.1 






^0 nnx 




146 


64 


9.6 


3.6 


149 


0.5 


47 


8.7 


58.936 


0.675 


147 


65 


9.1 


4,« 


9 


0,? 


78 


0 


41.154 


0,^*1 


148 


54 


9,4 


4.5 


1$ 


0.2 


9§ 


5.7 


69.583 


0.683 


149 


4? 


9.2 


4.1 


J? 


0.3 


143 


8.7 


34.895 


0.70? 


150 


75 


10 


?,9 


10 


0.2 


?S 


0 


39,545 


o.w 


ISl 


57 


9.4 


3.9 


10 


0.2 


54 


3.4 


51.111 


0.$92 


1 m 


47 


9.4 


4 


42 


0.3 


$5 


8.6 


51,»4$ 


0.739 


1 i» 


70 




3.5 


$ 


0.1 


« 


5 


37.231 


0.621 


154 


59 


9.2 


4 


9 


0.2 


49 


4.4 


47.143 


0.$90 


155 


W 


9.2 


4.1 


12 


0.1 


56 


2.8 


60.357 


0.680 


156 


49 


9.1 


3.6 


59 


0.2 


35 


0 


54.571 


o,7« 


157 


59 


9.5 


4.1 


7 


0.1 


92 


6.5 


45.978 


0.660 


158 


66 


9,2 


?,? 


1 


0,2 


46 


0 


49.OOQ 


0.660 


159 


57 


9.2 


4.3 


1^ 


02 


?1 


5.4 


53.704 


0.679 


\W 


41 




4.6 


306 


1.2 


46 


0 


75.217 


0.763 


161 


3? 






4 


0.2 




5.1 


69.265 


0.758 


162 


35 








0.4 




3.3 1 64.561 





SUBSTITUTE SHEET (RULE 26) 



wo 96/12187 9 / 2A FCTAJS95/01379 



Testing Data for Osteoporosis 
Diagnostic System 



1 PATrerrr 




CAL 


PHOS 


ETWO 


PROQ 


TALP 


TINTES . 


PLFYER 


SCORE 


1 123 




9.4 


4.2 


-1 


0.3 


61 


0 


46 


0.722 


1 128 


27 


9A 


4.4 


■I 


0.4 


58 


0 


43 


0.862 


131 


39 


9.1 


3.9 


156 


11.5 


44 


0 


59 


0.820 


173 


26 


9.5 


3.8 


41 


0.2 


98 


0 


60 


0.790 


201 


30 


9.1 


3.5 


123 


3.3 


90 


0 


53 




206 


39 


9.2 


2.6 


5Q 


1.5 


74 


89 


42 


0.729 


213 


46 


9.7 


3.4 


37 


0.2 


90 


0 


36 


0.708 


220 


46 


8.8 


3.4 


118 


4.3 


90 


9 


63 




223 


49 


?.7 


3,3 


72 


0.4 


78 


8.6 


38 


0.717 


234 


38 


9.5 


3.2 


186 


12,$ 


50 


0 


46 


0.803 


243 


2? 


9.3 


4.4 


65 


0.3 


58 


0 


40 


0.857 


244 


5? 




3,7 


9$ 


Q,3 


58 


0 


49 


0.725 


245 


30 


?.5 


3.? 


-1 


0.1 


64 


0 


53 


0.800 


247 


28 


9.1 


3.6 


■\ 


0.3 


58 


0 


50 


0.821 


249 


22 


9.3 


3.5 


42 


0.2 


98 




38 


0.782 


253 


24 


9.2 


3.4 


•1 


Q.3 


$8 


0 


63 


0.797 


1 254 


30 


9.6 


4 


H3 


10.1 


45 


0 


31 


0.894 


1 256 


19 


9.1 


3.8 


219 


Or? 


53 


6.1 


32 


0.871 


II 258 


30 


9.3 


3.4 


60 


0.6 


68 


0 


33 


9.929 


II 259 


21 


9.2 


4.4 


44 


0,2 


108 


0 


54 


0.805 


260 


23 


9.1 


3.8 


-1 


0.3 


62 


0 


61 


Qrm 


261 


23 


9.4 


3.4 


-1 


0.2 


85 


0 


32 


0.795 


265 


24 


9.4 


3.2 


191 


0.2 


71 


0 


30 


9,829 


266 


2<} 


9 


3.9 


91 


0.3 


51 


5,5 


47 


o,?<59 


268 


20 


9.3 


4.3 


62 


0.2 


61 


0 


60 


0.857 


269 


19 


8.9 


3.9 


53 


0,8 


68 


0 


42 


0.837 


271 


20 


8.9 


4.1 


277 


18.9 


65 


0 


40 




272 


24 


9 


3.3. 


58 


A A 




n 


dO 


0 815 


274 


19 


9.7 


3.8 


44 


0.3 


120 


0 


70 


0.776 


276 


2S 


9.2 


3.9 


-1 


0.9 


65 


5 


27 


0.837 


279 


21 


IQ 


4.6 


115 


0.5 


106 


0 


70 


0,834 


280 


21 


9.5 


4.1 


292 


0.9 


68 




S7 


0.856 




20 


9.1 


3.^ 


-1 


0.3 


47 


0 


53 


0.831 


284 


27 


9 


3.3 


37$ 


0,$ 


43 


0 


37 


0.855 


285 


31 


9.6 


3.7 


208 


1.2 


60 


0 


62 


0.842 


286 


27 


9 


3.2 


350 


1.5 


65 


0 


60 


0.819 


287 


22 


9 


3.1 


45 


0,2 


66 


0 


39 


9 805 


290 


23 


. 9.2 


4.3 


57 


0.3 


104 


0 


47 


0,817 


292 


24 


9.4 


3.7 


35 


0.3 


62 


0 


43 


0.?29 1 


296 


21 


9.2. 


3.? 


65 


0.3 


111 


0 


55 


0.767 


298 


26 


9.4 


3.2 


248 


0.8 


57 


0 


35 


0.842 


321 


2? 


8.9 




84 


0.2 


66 


0 


42 


0,7?5_ 




22 


9.4 


3.6 


45 


0.4 


50 


0 


37 


0.847 


331 


48 


9.7 


3,1 


47 


0.3 


84 


0 


58 


0.677 


332 


51 


9.8 


3.5 


-I 


0.3 


64 


0 


40 


0.707 


333 


5i 


9.2 


3.3 


166 


1.2 


63 


0 


56 


0.705 


340 


^ 




U. 






SL 


0 


rJd. 


0 728 



SUBSTITUTE SHEET (RULE26) 



wo 96/12187 



1 0 / 2 A PCT/US95AH379 



PATIENT 


AGE 




PH9S 




PROG 


TALP 


TINTES 


PUYER 






50 


9.6 


3.8 


76 


0.3 


$2 


0 


71 


0.730 


346 


44 


9-2 


3.4 


111 


4,9 


72 


12.4 


. 59 


0.977 


347 


41 


9.5 


3.7 




0.2 


77 


0 


69 


0.735 


359 


54 


9.3 


3.6 


95 


0.2 


77 


0 


. .. $?■ 


0-993 


?W 


34 


9.4 


3.6 


$2 


0.4 


76 


0 


.. $3 


. ,0.721. 


364 




HQ 


4,2 


104 


0.3 


62 


9 


57,. 


0.8U 


365 


46 


? 


2r5 


97 


6.7 


54 


0 


64 


0.698 


375 


Si 


?.3 


3,9- 


75 


„ . 0.2 


62 


0 




0.718 


382 


4? 


8.7 


4.8 


36 


0.2 


63 


0 


66 


0.796 


386 


24 


9.3 


3.8 


-1 


0.3 


«8 


0 


65 


0.812 


389 


22 


9.1 


4.3 


-1 


0.3 


« 


0 


« 


0.831 


391 


32 


9.7 


3.3 


-1 


0.4 


57 


0 


§4 


0.797 


393 


29 


9.6 


3.5 


30 


0.5 


47 


0 


$1 


0.828 


?94 


27 


8.9 


4.1 


-1 


0.4 


45 


0 


5$ 


0.848 


402 


38 


8.7 


3.5 


33 


0.6 


43 




46 


0.786 


403 


?2 


9.2 


3.4 


-1 


0.3 


21 


0 


...21. 


0,842 


404 


49 


9.2 


4.5 


-1 


0.1 


100 


0 


53 


0.73$ 


407 


52 


9,4 


33 


152 


0,4 


71 


0 


52 


0.704 


408 


24 


9 


3,7 


^32 


0.$ 


87 


0 


49 


P 819 


409 


46 


9.2 


4.1 


69 


7.1 


52 


0 


58 


0.698 


412 


52 


9.7 


3.§ 


-1 


0.2 


$8 


0 


29 


0.704 


416 


44 


9.7 


3.4 


142 


0.7 


80 


0 


79 


0.692 


417 


49 


9.5 


3.4 


-I 


0.2 


75 


0 


47 


0.708 


418 


49 


\Q 


4.2 


403 


1.3 


88 


0 


75 


0.678 


419 


2i 


?.3 


3 


-1 


0.4 


82 


0 


65 


0.767 


422 


38 


9.2 


3.3 


198 


0.2 


76 


0 


?4 


0.780 


539 


4$ 


?.3 


4.5 


793 


0.7 


46 


0 


33 


0.813 


545 


26 


9.7 




96 


0.3 


65 


0 


42 


0.825 


549 


26 


9.5 


3.9 


50 


0.4 


46 


7.2 


49 


0.8^2 


560 


32 


8.? 


3.4 




0.2 


79 


7.4 


$1 


0.791 


562 


23 


9 4 


3 2 


-1 


0 4 


58 


0 


36 


0.815 


579 


20 


9 


3.4 


4? 


0.5 


80 




53 


0.797 


588 


3$ 


?.3 


3.5 


-1 


0.1 


4? 


9 


51 


0.793 


591 


28 


9.7 


3.7 




0.5 


?2 


0 


54 


0.838 


593 


35 


9 


3.7 


108 


0.3 


$? 




41 


0.811 


597 


32 


9.3 


3.5 


-1 


0.3 


53 


0 


71 


0.801 


602 


49 


9.3 


3.6 


180 


4.3 


$9 


Q 


46 




2?20 


29 


?.2 


4.4 


?2 


0.2 


49 


0 


41 


Q,S77 


3000 


28 


9A 


3.7 


49 


0.4 




0 


53 


Q,838 


3020 


25 


9.6 


4 


164 


16.8 


56 


0 


41 


Q.m 


100 


52 


8.7 


3.1 


371 


0.3 


59 


0 


54 


9793 


113 


53 


9.8 


3.5 


3 


0.2 


$3 


0 


$7 


0.690 


12$ 


42 


9 


4 


99 


l$,4 


83 


14.4 


26 


0984 


130 


50 


?.3 


3.5 


33$ 


o.e 


49 


0 


25 


9.754 


139 


74 


8.9 


3.7 


?09 


0.4 


35 


5.9 


36 


9fr48 


160 


60 




3.9 


SI 


0.3 


47 


5.2 




0 -^04 


176 








Hi. 


QJ_ 




,. g 


a_ 


0 h77 




SUBSTITUTE SHEET (RULE 26) 



wo 96/12187 4 « / PCTAJS95/01379 



11/24 



PATIENT 


AGE 


CAL 


PHOS 


ETWO 


PRTC 


TALP 


TINTES 


PLIVER 


SCORE 


177 


39 


9A 


3.1 


71 


0.2 


82 


0 


45 


0.751 




42 


9.1 


2.4 


177 


15.6 


54 


0 


70 


0.757 


211 


50 


9,5 


3,1 


69 


0.2 


77 


0 


60 


0.696 


214 


25 


9.2 


3.7 


-1 


0,9 


43 


0 


29 


0.852 


251 


24 


9.5 


3.9 


-1 


0.1 


61 


0 


25 


0.844 


273 


?<? 


9.4 


3 


99 


0.2 


156 


8.3 


26 


9,77§ 


299 


27 


9.2 


3.2 


41 


0.2 


77 


0 


58 


0.788 


m 


?5 


9 


3.9 


55 


0.3 


90 


8.7 


62 


0.773 


316 


41 


9.4 


3.5 


75 


7,3 


44 


0 


27 


0.733 


320 


52 


9.3 


3,4 


57 


0.2 


58 


0 


55 




355 


54 


9.1 


3.6 


93 


0,4 


55 


0 


40 


0.718 


358 


54 


9.5 


3.6 


38 


0.3 


71 


0 


54 


0.688 


361 


49 


9 


3.8 


31 


0.3 


74 


9.5 


50 


0.744 


?71 


49 


9,5 


3,8 


-1 


0.2 


72 


10.3 


42 


(>m 


377 


52 


9.6 


2<J 


120 


0.2 


69 


0 


53 


0.717 


385 


55 


9 


3.3 


-1 


0,2 


97 


0 


73 


0.646 


398 


53 


9,3 


4.6 


30 


0.2 


97 


10.3 


28 




406 


26 


9.6 


3.4 


25 


0.3 


62 


0 


79 


0.806 


414 


52 


9.3 


3.9 


56 


0,4 


66 


P 


72 


0.710 


421 


21 


9 


3.5 


72 


0,3 


69 


0 


52 


0.818 


536 


48 


9.4 


3.4 


54 


3.2 


56 


0 


63 


om 


544 


38 


9.3 


2.4 


126 


0.2 


64 


0 


59 


0.747 


566 


36 


9.3 


4 


-i 


0.5 


41 


0 


44 


Q.m 


603 


27 


8.9 


4.4 


42 


0.3 


48 


0 


38 


0.872 


109 


54 


9.4 


3.3 


64 


0,3 


121 


9.6 


65 


0.652 


110 


58 


9.4 


3.9 


20 


0.1 


85 


0 


58 


0.663 


129 


45 


8.8 


3.7 


197 


10.8 


56 


0 


53 


0.733 


141 


4$ 


9,7 


3,6 


HJ 


2,1 


68 


0 


44 


0.709 


158 


45 


8.9 


3.3 


243 


0.8 


63 


0 


43 


0.763 


164 


51 


9.3. 


3,5 


139 


0,5 


59 


0 


43 


0.727 


1 /U 




Q A 






A 9 




n 




0 7'^7 


178 


36 


9.7 


3.2 


139 


0,4 


40 


0 


46 


0.824 


180 


38 


8.7 


3.4 


48 


0.4 


52 


0 


39 


0.783 


m 


?4 


8.7 


3.4 


78 


11.1 


67 


6.8 


3J 


0.§17 


183 


44 


9. 


3.6 


131 


0.2 


35 


0 


54 


0.807 


\9l 


45 


9.3 


3.5 


-i 


0.4 


45 


0 


58 


0.760 


194 


29 


9.2 


3.7 


89 


0.1 


82 


0 


5? 


0.824 


200 


45 


9.2 


3.7 


226 


1.6 


46 


0 


27 


0.786 


204 


25 


9.5 


3.5 


170 


0.3 


43 


0 


46 


Q,m 


29? 


42 


9^ 


2.9 


539 


0,4 


50 


0 


64 


0-7J9 


216 


4$ 


9.2 


3.6 


(56 


0.5 


43 


0 


39 


0.780 


218 


43 


9.5 


4.4 


133 


0.7 


59 


0 


41 


0.791 


225 


J7 


9 


2.7 


159 


1.4 


64 


6.7 


48 


0.766 


226 


72 


9.9 


4.9 


46 


0,3 


77 


8.9 






227 


4^ 


9.6 


4.4 


-1 


0.2 


94 


0 


38 


0.757 


230 


43 


9.3. 


3.4 


94 


0 


71 


0 


77 


0.721 


22L 






3.2 


-» 


0.3 




SL 




0.777 



SUBSTITUTE SHEET (RULf 26) 



wo 96/12187 1 2 /2 4 PCT/US95/01379 



PATiErn- 


AGE 




PH<?5 






TALP 


TINTK 


PLIVER 


^ 

SCORE fl 


233 


50 


9.1 


3.5 


43 


0.3 


73 


0 


23 


0.709 


238 


4^ 


8.3 


3,4 


234 


0.2 


55 


5.9 


46 


1.094 


239 


40 


9.5 


3.§ 


43 




71 


0 


63 


0.758 


246 




9,8 


3.3 


-1 


0.3 


42 


0 


58 


0.804 


248 


51 


8.9 


3.8 


76 


0.2 


$1 


0 


62 


0.724 


255 


19 


9.$ 


4.3 


49 


0.3 


86 


0 


34 


0.841 


257 


23 


8.9 


2.9 


118 


2.7 


75 


0 


21 


0,&1^ 


262 




9.1 


3.$ 


-I 


0.2 


47 


0 


24 


0.844 


267 


?3 


9 


3.8 


170 


0.4 


84 


0 


50 


0.825 


1 ?70 


21 


9.1 


3.9 


47 


0.9 


40 


0 


70 


0.850 


1 277 


27 


9,4 


4.§ 


58 


0.2 


138 


0 


69 


0.811 


1 278 


29 


9.4 


3.9 


404 


1.2 


« 


7.9 


46 


0.864 


282 


28 


9-4 


3.9 




0.2 


57 


0 


36 


0.859 


283 


25 


IQ 




-1 


0.3 


48 


0 


32 


0.847 


288 


24 


9.7 


4.1 


65 


0.3 


107 


7,2 


35 


0.816 


289 


27 


9.4 


3.4 


$8 


0.6 


87 


0 


37 


0-812 


291 




9.$ 


3.3 


-1 


0.3 


58 


0 


28 


0.823 


293 


26 


9.6 


4 


32 


0.2 


68 


0 


35 


0.838 


295 




9.1 


3.3 


57 


0.5 


51 


0 


63 


0.324 


310 


47 


9.5 


3.6 


306 


1.5 


53 


0 


46 


0-730 


314 


48 


9.1 


2.7 


-1 


0.2 


64 


0 


53 


0.704 


328 


37 


9.4 


3.2 


54 


0,4 


74 


0 


55 


0.762 


335 


47 


9.2 


4.2 


49 


0,4 


74 


0 


56 


0.745 


341 


37 


9.3 


3.3 


-1 


0.1 


68 


0 


56 


0.758 


349 


50 


8.8 


3.1 


151 


0.1 


67 




20 


0.710 


353 


50 


9 


4.1 


374 


0.2 


53 


0 


30 


0.759 


.354 


58 


9.7 


4.3 


74 


0.3 


48 


6.8 


45 


0.726 


363 


50 


9.1 


3.3 


7$ 


0.$ 


93 


13.8 


31 


0.694 


366 


24 


9.4 


2.7 


-1 


0.5 


50 


0 


31 


0.806 


372 


21 


10 


2,9 


56 


0.7 


79 


0 


61 


0.789 


378 


42 


9.2 


3.5 


169 


A 1 

0.1 


IJL 


t\ 
V 




n 


381 


62 


9.5 


3.$ 


56 


0.1 


63 


0 


30 


0.672 


384 


« 


9.4 


4 


34 


0.2 


98 


25.3 


33 


0.644 


387 




9A 


3.7 


38 


0.2 


133 


0 


55 


0.669 


388 


$4 


9.1 


3.1 


87 


0.3 


80 


0 


75 


0.636 


396 


57 


8.8 


3.1 


7$ 


0,2 


83 


7.7 


51 


0.659 


397 


58 


9.§ 


3.3 


73 


0.3 


85 


0 


52 


0,665 


400 


50 


9.1 


3.5 


-1 


0.2 


99 


9.7 


50 


0.671 


410 


37 


9.5 


2.8 


-1 


0.3 


44 


0 


48 


0.769 


411 


49 


9.4 


3.1 


138 


7.5 


69 


0 


27 


0.674 


415 


, 3A, 


9.« 


4.6 


-I 


0.4 


42 


0 


72 


0.861 


424 


54 


9.4 


3.1 


59 


0.4 


HOI 


0 


35 


0 662 


507 


45 


9 


3.4 


46 


0.3 


55 


0 


65 


0.745 


513 


25 


9.4 


3.5 


52 


0.2 


72 


0 


52 


0,815 


533 


61 


9.6 


3.9 


-1 


0.1 


80 


0 


61 


0.655 


535 


24 


9 


3.3 


-1 


0.3 


68 


0 


55 


P,795 




48 




3.2 


32 


QJ_ 


— 


0 


58 


0.679 



SUBSTITUTE SHEET (RULE 26) 



wo 96/12187 



13 



PCTAJS95/01379 



PATIENT 


AGE 


CAL 


PHOS 


ETWO 


PROG 


TALP 


TINTES 


PLIVER 


SCORE 


540 


45 


9.5 


2.9 


391 


0.7 


49 


7 


35 


0.955 


541 


54 


9.3 


4.4 


-1 


0.2 




0 


37 


9,7|4 


542 


32 


9.2 


3.4 


80 


0.6 


51 


0 


60 


9,8?9 


546 


62 


9.3 


4.5 


100 


0.1 


45 


9 


59 


0.711 


548 


52 


9.2 


4.1 


-1 


0.4 


55 


0 


59 


9.719 


550 


48 


9.9 


3.7 


92 


0.3 


127 


0 


78 




554 


36 


9.1 


3.5 


42 


0.4 


44 


9 


37 


0.808 


561 


30 


8.9 


3 


54 


0.2 


55 


0 


54 


P.899 


568 


53 


9.6 


3.2 


-1 


0.3 


49 


0 


40 


9.599 


569 


26 


9.2 


4.5 


274 


0.5 


80 


9 


30 


0.870 


571 


45 


9.1 


3.3 


402 


1 


72 


0 


54 


0.731 


572 


58 


9.5 


3.9 


52 


0.2 


57 


7.3 


23 


0.700 


576 


65 


9.5 


3.3 


43 


0.2 


78 


9 


56 


0.632 


580 


50 


9.3 


3.2 


60 


0.2 


76 




43 


(>m 


582 


44 


8.9 


3 


119 


0.2 


80 


7.1 


24 


9,9§7 


583 


45 


8.9 


3.2 


353 


0.2 


68 


0 


64 


0.744 


586 


41 


8.8 


3.7 


223 


0.7 


73 


12.4 


22 


1.206 


589 


52 


9.3 


3.5 


-1 


0.3 


88 


0 


45 


0.678 


592 


55 


9.2 


3.7 


95 


0.3 


50 


0 


57 


0.717 


595 


46 


9.1 


3 


128 


4.5 


90 


8.2 


35 


1.000 


596 


46 


9.5 


3.6 


-1 


0.2 


53 


0 


47 


0.757 


604 


56 


9.3 


4.3 


-1 


0.2 


84 


0 


53 


9,$?4 


605 


47 


8.9 


3.4 


-1 


1.4 


64 


0 


40 


0.740 


608 


46 


9.7 


3.8 


-1 


0.1 


69 


0 


48 


0.742 


2500 


30 


9.2 


3.6 


70 


0.2 


62 


19.1 


42 


0.842 


3010 


24 


9.4 


3.8 


178 


15.9 


33 


0 


34 


P.894 


106 


52 


9.1 


3.4 


44 


0.2 


86 


9 


55 


0.674 


117 


61 


9.5 


3.2 


-1 


0.1 


85 


0 


60 


0.634 


121 


46 


9.5 


2.9 


48 


0.3 


80 


0 


54 


0.690 


133 


47 


9 


4.2 


203 


15 


57 


7 


52 


0.939 


151 


36 


9.1 


3.8 


-1 


0.3 


80 


0 


57 


0.768 


156 


24. 


9.$ 


3.9 


114 


17.9 


41 


0 


37 




166 


SI 


?.2 


3.8 


60 


0.1 


60 


0 


42 


0.727 


167 


49 


8.5 


3.8 


7$ 


0,5 


97 


^9 


48 


9.?19 


179 


32. 


?.5 


3.^ 


59 


0.2 


61 


0 


47 


0.826 


182 


..Si. 


8.4 


3.7 


99 


9,2 


«4 


9 


31 


0.709 


182- 


Si 


9.2 


4.2 


-1 


9,1 


68 


8 


47 


0.705 


212 


37 


9.4 


3.7 


162 


18.6 


55 


0 


57 


9,«18 


217 


«7 


9.4 


3.2 


-1 


9,2 


«4 


6.3 


33 


0.628 


219 


S<L 


9.6 


3.7 


-1 


0.3 


81 




18 


0.707 


222 


Si 


9.1 


3.9 


30 


9.2 


102 


0 


61 


0,«7 


228 


SI 


9.4 


3.5 


141 


9.4 


$1 


0 


28 


9.739 


22i. 


41 


9 


4.3 


96 


0.3 


45 


0 


39 


0.868 


252 


52 


9.5 


3.3 


34 


0.2 


79 


0 


48 


0.681 


264 


28 


9.1 


3.4 


-1 


9.? 


50 


0 


«9 


9.5I« 


275 


25 


9.1 


i.s 


-1 


9,4 


38 


0 


34 


0.806 


1 297 






4.4 


177 


12^ 


72 


(L 


55 


0.873 



M 7 

SUBSTITUTE SHEET (RULE 26) 
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PATIENT 


ACE 


CAL 


nios 


ETWO 


PROG 


TALP 


TINTES 


PLIVER 


S?OB£.. 


317 


47 


9.$ 


2.4 


-1 


SL2- 


. $0 


0 


27 


0.709 


324 


4$ 


9.4 


4 


1 17 


0.4, 


75 


0 


54 


0.73? 


327 


4$ 


9 


3.4 


322 


1 


46 


0 


52 




329 


49 


9.5 


3.7 


n 


2.4 


78 


0 


45 




?37 


4$ 


9.3 


3.5 


m 


19.3 


48 


0 


42 


0.700 


352 


56 


8.8 


3.6 


73 


0.2 


53 


0 


35 


0.707 


357 


35 


9.4 


2.9 


89 


2.3 


71 


7,$ 


23 


0.79J 


367 


42 


9.1 


3.2 


58 


0.5 


90 


9.8 


64 


0.846 


369 


$? 


9.$ 


3.6 


31 


0.2 


8$ 


0 


39 


O.W 


379 


42 


8.$ 


3.4 


156 


$.3 


84 


0 


61 


0-742 


380 


49 


9.6 


4.2 


-l 


0.2 


71 


u 


52 




383 


62 


9 


3.5 


81 


0.1 


66 


7 


52 


0.665 


392 


63 


9.5 


3.4 


-1 


0.3 


98 


0 


76 


0.620 


517 


25 


10 


4 


31 


0,1 


50 


0 


50 


0.853 


5?$ 


47 


9.1 


3.8 


94 


0,2 


$2 


0 


51 


0.750 


543 


59 


9.4 


4.2 


-1 


0.2 


104 


0 


20 


0.662 


553 


48 


8.8 


4,1 


35 


0.5 


57 


0 


66 


0.759 


556 


69 


9.6 


4 


77 


0,1 


44 


6.8 


53 


0.668 


5§7 


29 


9.4 


4.? 


-1 


0.3 


67 


9.7 


33 


0.841 


575 


45 


9.1 


3.1 


406 


0,$ 


104 


0 


54 


0.704 


594 


57 


9.? 


3.2 




0.2 


79 


0 


40 


0.674 


599 


49 


8.9 


4 


42 


0,2 


79 


9 


51 


0.748 


2510 


2$ 


9.5 


3.9 




1$.2 


55 


0 


25 


Q.881 


1Q5 


75 


9.4 


4.1 


42 


0.4 


88 


8 


52 


0.605 


11$ 


$7 


8.9 


4 


-1 


0.3 


5$ 


0 


70 


0.647 


118 


55 


8,9 


4.3 


38 


0.2 


78 




37 


0.699 


119 


56 


9.3 


4,3 


75 


0.3 


41 


0 


43 


0.740 


127 


73 


9 


4.J 


30 


0.4 


54 


$.5 


35 


0.632 


134 


50 


9.4 


3.7 


486 


0.4 


54 


0 


$3 


0.737 


13$ 


3$ 


9.1 


3.$ 


J$7 


5,7 


154 


0 


74 


0.758 


140 


56 


0 




1 


M.i 




C\ 

V 




0 ftf%'\ 


H3 


55 


9.1 


4 


-\ 


0,2 


$5 


0 


38 


0.698 


145 


49 


9.3 


3.3 


67 


8,9 


« 


0 


54 


0.656 


147 




.2.L. 


3.$ 


75 


0.2 


103 


4.9 


78 


0.$70 


157 


52 


9.2 


3.5 


65 


0.2 


70 


0 


62 


0.701 


15? 


57 


9.4 


3 


-1 


0.2 


129 


30.2 


27 


o.wi 


161 


74 


9.3 


3.9 


-1 


0.1 


62 


0 


34 


0.623 


162 


62 


9 


4.1 


81 


o.> 


54 


7.2 


21 


0.698 


i» 


46 


9.2 


?,9 


128 


5.7 


$3 


0 


35 


0.704 


1$5 


57 


9.2 


2.9 


-1 


0.2 


108 


0 


39 


0.629 


1«? 


59 


9.5 


3.9 


$1 


0,1 


$4 


0 


34 


O.W 


174 


40 


9.5 


3.2 


261 


4.6 


85 


0 


22 


0,7^ 


175 


39 


9.1 


2.3 


m 


3.5 


90 


0 


36 


0.726 


185 


78 


8.7 


2.9 


-1 


0.1 


99 


8.5 


35 


05« 


18$ 


73 


9.2 


2.9 


$2 


0.2 


55 


0 


42 


0.616 


187 


59 


9 


3.8 


55 


1.8 


70 


0 


37 


0.682 


19Q_ 






^ 
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SUBSTITUTE SHEET (RULE 26) 
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, PATIENT 


ACE 




PHOS 


ETWO 




TALP 


TtNTES 


PLJVER 


SCORE 


192 


42 


9.5 


4.1 


4« 


9.4 


58 


0 




0.767 


195 


59 


9 


3,4 


236 


0.1 


61 


0 


40 


0.688 


202 


67 


9,4 


3.2 


66 


0.3 


59 


0 


61 


0.641 


210 


«9 




3,2 


82 


0.1 


49 


4.7 


5$ 


9r$8l 


215 


45 


9.2 


3,8 


135 


1.8 


70 


0 


56 


0.729 


221 


52 


9.5 


4 


83 


9,8 




0 


57 


0.735 


224 


53 


8.7 


3.6 


110 


9,2 


52 


0 


27 


0.726 


229 


51 


9.5 


4.1 


33 


9.3 


97 


9 


51 


9.693 


237 


5? 


9.4 


4 


-1 


0.2 


90 


9.9 


34 


0.660 




48 


9.6 


3.9 


44 


9.5 


$8 


9 


31 


0.747 


263 


23 


9.6 


4.2 


39 


9.5 


«9 


0 


28 


0.853 


294 


25 


9.3. 


3.8 


91 


2.9 






54 


0.853 


308 


54 


8.8 


3.1 


-I 


0.2 


40 


0 


48 


0.692 


313 


53 


9.8 


4.3 


-1 


0.4 


90 


0 


40 


0.700 


319 


31 


9 


3.5 


296 


2.8 


49 


0 


27 


0.859 


322 


42 


9.4 


3,3 


-\ 


9,3 


$3 


^,6 


41 


0.757 


?2? 


5$ 


8.8 


3.6 


38 


0.1 


54 


6.5 


54 


0.685 


334 


$8 


8.8 


4,1 


83 


9.4 


44 


0 


39 


9,$7§ 


338 


42 


9,9 


3,2 


35 


I 


92 


9.4 


41 


0.719 


339 


52 


9,9 


3,7 


?4 


9,2 


78 


0 


« 


0.7U 


J4J 


44 


9.1 


3.4 




U 


75 


8 


55 


0.834 


344 


61 


9.2 


4.3 


82 


9,3 


116 


8,3 


65 


0.655 


m 


46 


9,2 


3.9 


140 


0.3 


58 




51 


0.905 


362 


53 


9.6 


3.9 


94 


0.1 


83 


0 


72 


0.703 


373 


« 


9.9 


3,4 




9,2 


83 


9 


66 


0,$54 


376 


70 


9.9 


4.7 


-1 


0.1 


109 


0 


38 


9<?27 


. 390 


72 


9 


3,8 


-I 


9,2 


198 


14.5 


28 


0.597 


395 


54 


9.2 


4.1 


-1 


0.3 


9$ 


0 


67 


9-$74 


399 


49 


9.4 




-1 


0.4 


64 


9 


59 


9,791 


401 


83 


8.7 


3.8 


31 


9.2 


111 


0 


55 


0.558 


413 


46 


9.4 


3.$ 


99 


9.1 


98 


25.1 


54 


0.746 


423 


44 


9.2 


3.6 


62 


0.2 


47 


0 


45 


0.7?5 


527 


3? 


9.3. 


3,3 


40 


9.4 




4.1 


5$ 


0.763 


552 


54 


9,4 


3,7 


-1 


0.2 


80 


0 


57 


0.679 




SO 


9,5. 


3 


108 


9.2 


12$ 


7.4 


50 




557 


59 


9.4 


3.5 


192 


0.2 


42 


9 


^ 


0.697 


5^3 


44 


9.5 


3.8 


37 


0.5 


4§ 


0 


39 


0.790 


565 


70 


9 


3,2 


-1 


9.2 


54 


9 


45 


0.623 


570 


?4 


9.2 


3,5 


59 


9.5 


53 


7 


32 


0.842 


573 


$5 


9.9 


4.2 


-1 


9.3 


59 


0 


56 


0.665 


574 


« 


9,2 


3.2 


122 


9,3 




0 


58 


0.649 


577 


46 


9.5 


4.1 


397 


.04 


56 


0 


70 


0.738 


585 


« 


9 


3.5 


56 


9.2 


77 


9 


21 


0.658 


590 


65 


9.2 


3.7 


-1 


9.2 


72 


0 


57 


0.640 


594 


55 


9.9 


3,$ 


-1 


9.2 


43 


9 


52 


0.705 


601 


45 


8.9 


3.1 


81 


3.4 


75 


6.6 


54 


9.??? 


607 




_2Jn 




-1 






0 








SUBSTITUTE SHEET (RULE 26) 



wo 96/12187 



PCTAJS95/01379 



PATIENT 


AGE 


CAL 


PHOS 


ETWO 




TALP 


TINTES 


PLIVER 




III 




9 


3.3 


34 


0.4 


68 


9.3 


55 


QMS 


124 


51 


9.2 


3.7 


35 


0.2 


$9 


0 


48 


0.704 


149 


52 


9,2 


3.9 


-I 


0,3 


82 


0 


72 


0.§8$ 


152 


56 


9 


2.5 


79 


0.2 


77 


0 


43 


0.659 


168 


5$ 


9.1 


4.1 


-I 


O.I 


76 


0 


45 


0.686 


196 


$1 


9.3 


3.3 


69 


0.1 


97 


8.8 


36 


0.642 


330 


$3 


8.9 


3,8 


-1 


0,1 




9.$ 


54 


0.648 


345 


65 


8.9 


4.2 


40 


0.1 


50 


5.3 


19 


0,§7$ 


356 


SI 


9.7 


3.5 


54 


0.2 


99 


0 


$0 


0,$73 


534 


66 


8.9 


3.5 


78 


0.1 


$1 


0 


53 


0.654 


SSI 


67 


9.3 


3,5 


47 


0.2 




10 


47 


0,M5 


SS9 


38 


9.3 


3,3 


137 


7,9 


85 


0 


66 


0.765 




51 


8.8 


3,9 


34 


0-3 


50 


0 


54 


0.722 




53 


9.3 


3.4 


62 


0.2 


5§ 


0 


54 


0.707 


104 


63 


9.3 


4.3 


77 


0.3 


88 


7.6 


41 


0.669 


UQ 




9,3 


3,8 


35 


0,2 


81 


0 


$1 


0.633 


125 


$3 


9.3 


4 


-1 


0.4 


92 


0 


55 


0.642 


138 


60 


9,3 


2,8 


58 


0.3 


84 


0 


77 


0,$34 


144 


68 


9.5 


4.6 


-1 


0.4 


118 


0 


30 


0.632 




$7 


9.8 


3.3 


-1 


0.2 


74 


0 


59 


0.623 


m 


44 


9.$ 


3.8 


65 


0.2 


122 


0 


44 


0.716 


m 


50 


10 


3.2 


-1 


0.3 


§7 


0 


24 


0.705 


\99 


36 


9.4 


3.2 


262 


0.6 


101 


14.1 


18.3 


0.761 


TXfl 


57 


9.7 


3.5 


33 


0.3 


74 


0 


57 




232 


§5 


9,4 


4 


80 


0,2 


§9 


0 


53 




236 


48 


9 


3.6 


73 


8.9 


43 


0 


33 


0.704 


240 


53 


9.1 


4.1 


41 


0.2 


5§ 


5 


39 


0.721 


242 


57 


9.4 


3.8 


-1 


0.1 


7| 


9.8 


34 


0.677 


315 


64 


9.7 


3,5 


4$ 


0,2 


83 


11,4 


37 


0.639 


m 


49 


9.8 


3.2 


168 


3 


4« 


0 


34 


0.701 


325 


57 




4.3. 


_-l 


0J_ 


124 


0 


43 


Ow2 


336 


55 


9.2 


3.6 


-1 


0.3 


76 


10 


28 


0.676 


348 


58 


8J., 


3,5 


-I 


0.2 


55 


0 


55 


0.^72 




62 


9.4 


3.9 


-1 


0.2 


68 


0 


45 


0.663 




77 


9.5 


4.1 


-I 


0.4 


98 


10.2 


21 


0.596 


374 


, 7§ 


9.9 


3.7 


-J 


0.3 


78 


0 


54 


0.600 


528 


80 


9 


3,9 


6 


0.2 


103 


0 


71 


0,571 


529 




9.1 


3.5 


-1 


0,1 


§7 


0 


54 


0.634 


538 


79 


9,4 


3.9 


-1 


0.3 


9^ 


9,4 


39 


0.582 


547 


73 


9.3 


3.6 


126 


0.1 


§4 


0 


22 


o,«6 


559 


68 


9.7 


3.6 


91 


0.1 


37 


0 


57 


0,<?$9 


564 


80 


9.5 


3,$ 


33 


0.2 


113 


10.7 


52 


0.565 


598 


69 


9A 


4.4 


128 


0.7 


72 


2.5 


57 


0.660 




fr» 


9.2 


4,5 


-1 


0.4 


66 


0 


47 


0.671 


U 606 


78 


9.5 


3.9 


-1 


0,1 


77 


0 


58 


0,59$ 
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wo 96/12187 
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PCTAJS95/01379 



• • Program 
Author 
Date 
Purpose 



Calc. prg 

Terminal solutions (JJK) 

14 August 1992 

Quios calculating program 



set decimal to 16 
set fixed on 

****** Set CoefTicient infonnation to valiables ' 
0.000000000 to oumumb 
1009.33 to constant 



store 
store - 
store - 0.014786248 
store 424.462000000 
store - 66.3468 
store - 0.025761569 
store • 0.019715699 
store 0.027232306 
store - 0.000000263 
store . 12.0319 
store 7.499334732 
store • 2.220366330 
store 0.313008013 
store - 0.016909117 
store 0.001369629 
store - 0.000428831 
store 0.001798351 
store 0.001841052 
store 0.007649946 
store 0.002153604 
store 4.606320653 
store - 0.119876531 
store - 0.011786744 
store - 0.000170224 
store 0.OO0000706 
outnumb^ (constant) 
outnumb^outnumb + 
oumtmib=outnttmb 
oumiunb=outottmb 
oumiunb^oumumb -i- 
oumiimb»outnumb + 
oumiuib^oumiimb 
ouaumb»outmmib + 
outnumb=^outnurob + 
outnurob=outnumb + 
oumumb=outnumb + 
outnumb^utnumb + 
oumtunlm)utnimib + 
oumimib=outmmib + 
oumiunb^outnimib -*- 
outnumb=outnumb -** 
oumumb=oumumb 
outnumb=outnumb 
outnumb^outnumb 
outnumb=H>utnumb + 
outnumb^utnumb 
outnumb^umiunb + 
oumiunb^oumumb 
oumumb=outnumb + 



to c _age 
to c_ cal 
to c_cal2 
to c_csage 
to c_csphos2 
to c_cstintes 
to c_etwoimx2 
to c_ phos 
to c_ phos2 
to c_ phos3 
to c_ pbos4 
to c_ phos5 
to c_ prog2 
to c_ progpliv 
to c_ progtint 
to c_ tnetwoimx 
to c_ mtintes 
to c_ agephos 
to c_ caD 
to c_ caM 
to c_ talp 
to c_talp2 
to c.talp3 



(c_agc * age) 

(c_ cal * cal) 

(cl cal2 * (cal * cal)) 

(c_ csage * (cos(age))) 

(c_ csphos2 * (cos(phos * phos))) 

(c_ cstintes * (cos(tintes))) 

(c_ etwoimx2 * (etwoimx * etwoimx)) 



(c_ phos * 
(c_ phos2 * 
(c_ phos3 • 
(c_ phos4 • 
(c_ phos5 * 
(c_prog2 * 
(c_ progpliv * 
(c_ progtint 
(c_ metwoimx 
(c_ tntintes * 
(c_ agephos * 



phos) 
(phos * 
(phos * 
(phos ♦ 
(phos * 
(prog ♦ 
• (pr«>8 



phos)) 
phos. 
phos * 
phos ^ 
prog)) 
* pliv)) 



' phos)) 
phos * 
phos * 



phos)) 
l^os * 



phos)) 



(prog * tintes)) 
^ (tan (etwoimx))) 

(tan(tintes))) 
(age • phos)) 



(c_ caD 
(c_ cal4 
(c_ talp 
(c_ talp2 
(c_ talp3 



(cal • 
(cal 

talp) 
* (talp 
' (talp 



cal * 
cal * 

talp)) 
lalp * 



cal )) 
cal * 



talp)) 



cal)) 
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fi 4 \ = 0_N-1 4 input nodes 

Nl 8 j = 0_N1 I 8 neurons in first hidden layer 

N2 1 1 output neurons 

f(*) 1 Logistic neuron output function 

1+e 





100 




0 


Imax. = 


25 


Imin. = 


10 




40 




0 




40 




0 



Irange. = Imax. -Imin. 



Wl = READPRN (link l_wgt) 
W2 = READPRN (link 2_wgt) 



Wl = 



-0.26 


-0.33 


-0.15 


-0.44 




0.03 


0.47 


0.09 


O.QS 


-0,17 1 


0 


-1.16 


-0.07 


-0.12 


-0.05 


0.06 


1.25 


-0.26 


0.22 


-0.09 


0.01 


0.73 


0.21 


0.15 


0.02 


0.02 


0.74 


-0.1 


0.21 


0.07 




0.23 


1.18 


0.05 


-0.63 


0.18 


[a3i 


1.25 


0.27 


0.01 


0.07 



9x(5+l) weithi matrix between input 
layer and the first hidden layer 



W2 = (0.43 1.46 1.39 3.91 3.96 2.16 2.25 -6.36 -4.1) 2x9 weight matrix between 

hidden layer and output 
layer 
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